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ABSTRACT 

We present a comprehensive analysis of the whole sample of available XMM-Newton observations of High Mass X-ray Binaries 
(HMXBs) until August, 2013, focusing on the FeKcr emission line. This line is a key tool to better understand the physical properties 
of the material surrounding the X-ray source within a few stellar radii (the circumstellar medium). We have collected observations 
from 46 HMXBs, detecting FeKo- in 21 of them. We have used the standard classification of HMXBs to divide the sample in different 
groups. We find that: (1) different classes of HMXBs display different qualitative behaviour in the FeKff spectral region. This is 
specially visible in SGXBs (showing ubiquitous Fe fluorescence but not recombination Fe lines), and y Cass analogs (showing both 
fluorescent and recombination Fe lines). (2) FeKo- is centred at a mean value of 6.42 keV. Considering the instrumental and fits uncer¬ 
tainties, this value is compatible with ionization states lower than Fe XVIII. (3) The flux of the continuum is well correlated with the 
flux of the line, as expected. Eclipse observations show that the Fe fluorescence emission comes from an extended region surrounding 
the X-ray source. (4) We observe an inverse correlation between the X-ray luminosity and the equivalent width of FeKir (EW). This 
phenomenon is known as X-rays Baldwin effect. (5) FeKn is narrow ((T;,„(, < 0.15 keV), reflecting that the reprocessing material 
does not move at high speeds. We attempt to explain the broadness of the line in terms of three possible broadening phenomena: line 
blending, Compton scattering and Doppler shifts (with velocities of the reprocessing material V ~ 1000 km/s). (6) The equivalent hy¬ 
drogen column (Nh) directly correlates with the EW of FeKff, displaying clear similarities to numerical simulations. It highlights the 
strong link between the absorbing and the fluorescent matter. (7) The observed Nh in Supergiant X-ray Binaries (SGXBs) is in general 
higher than in Supergiant Fast X-ray Transients (SFXTs). We suggest two possible explanations: different orbital configurations, or 
a different interaction compact object - wind. (8) Finally, we analysed in more detail the sources ICR J16320-4751 and 4U 1700-37, 
covering several orbital phases. The observed variation of Nu between phases is compatible with the absorption produced by the 
wind of their optical companions. The obtained results clearly point to a very important contribution of the donor’s wind in the FeKcr 
emission and the absorption when the donor is a supergiant massive star. 
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1. Introduction 


Since the early stages of X-ray astronomy, Fe lines in the spec¬ 
tral region of ~6-7 keV (the Fe complex) have been studied in 
a large number of X-ray sources given its fruitfulness as a tool 
for plasma diagnostics. They were reported for the fi rst time in 
the supernova remnant Cas A (ISerlemitsos et al.ll973l) . and only 
two years later in a High Mass X-ray B inary (HMXB) using 
the Ariel 5 satellite (ISanford et al.lll97^ . The most recent X- 
ray space missions (Swift, Suzaku, Chandra and XMM-Newton) 
have triggered a notable improvement in the attainable spec¬ 
tral resolution and effective area, permitting to distinguish be¬ 
tween different emission features in the Fe complex: narrow and 


broad fluorescence lines (FeKtr and FeKjS), Com pton shoulders 
and rec ombination lines (Fe XXV and Fe XXVl) (ITorreion et al.l 
l2010bh . This improvement has given a remarkable impetus in the 
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study of the Fe complex, and motivates a comprehensive analy¬ 
sis in HMXBs. 

In particular, FeKo' has bee n proven as a fundamenta l 
tool in the study of HMXBs (iM artlnez- Nunez et ^ I20l4 
iRodes-Roca et al.ll201 lUvan der Meer et al.ll2005h . The origin of 
the fluorescence e mitting region h as been discussed by many 
authors in the past. iNagasel (Il989h considered accretion disks 
and the matter stagnated in the accretion and ionization wakes 
in the stellar wind a s plausible areas of FeKa production. 
IWatanabe et alJ (l2006h analysed the classical HMXB Vela X-1, 
and proposed the extended stellar wind, reflection off the stellar 
photosphere and an accretion wake, as the most likely candidates 
for fluorescence reprocessing regions. In any case, FeKa is very 
sensitive to the physical conditions of the vicinity of the X-ray 
source, and provides remarkable information that must be anal¬ 
ysed. 

Fluorescence is produced as a consequence of the X-ray il¬ 
lumination of matter. When an Fe atom absorbs a photon car¬ 
rying sufficient energy to remove an electron from its K-shell 
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(E>7.2 keV), the vacancy can be occupied by another electron 
from an outer shell. If the electron comes from the L-shell, the 
transition produces FeKa emission. FeK^S emission is produced 
when the vacancy is filled by a former M-shell electron. When 
Fe is more ionized than Fe XIX, t he fluorescence yield starts to 
decrease with the ionization state (iKallman et al.ll2()()4l) . There¬ 
fore, FeKtr is a footprint of not extremely ionized Fe (less than 
Fe XX). On the other hand, recombination lines Fe XXV and 
Fe XXVI unveil the presence of very hot gas, where Fe atoms 
are almost completely stripped. 

Previous comprehensive surveys of th e Fe c omplex in 
HMXBs were carried ou t bv iGottwald et al.l (Il995h using EX- 
OSAT and iTorreion et al.l (l2010bh using the High Energy Trans¬ 
mission Gratings (HETGS) on board Chandra. The high spec¬ 
tral resolution provided by Chandra gratings proved to be instru¬ 
mental in disentangling the different ionization species present 
in the Fe complex. However, the relatively low throughput of 
the instrument allowed to study only the brightest binaries. In 
this study we increase significantly the previous sample by us¬ 
ing the high throughput of XMM-Newton EPN. This has allowed 
us to include fainter systems (like Be X-ray binaries (BeXBs) or 
SFXTs in quiescence) while the moderate resolution of the EPN 
CCDs has allowed us to test previous correlations based on a 
small sample. 

HMXBs are specially susceptible to be studied using the 
Fe complex, on account of the significance of the circumstellar 
medium in the observable phenomena. These systems consists 
of a compact object, either a neutron star (NS) or a black hole 
(BH), accreting matter from a massive OB star (usually called 
optical or normal star of the system). In HMXBs the observed 
luminosity is commonly powered via accretion. Consequently, 
the way that matter is accreted from the donor directly defines 
the observable luminosity features of every source. 

When the optical star is a Be star, the system is a BeXB. Be 
stars are fast-rotating BIII-V stars, which have shown spectral 
emission lines at some point of their lives. They also show an 
excess of infrared emission, when they are compared to non Be 
stars of the same spectral type. These observables are explained 
appealing to an extended circumstellar decretion disk. BeXBs 
are usually transient in the X-rays, although some systems ex¬ 
hibit a persistent quiescence emission (F< 10^"^“^^ erg/s). The 
outbursts have been traditionally classified in two types. Type I 
outburst (F< 10^^ erg/s) are related to periastron passages. Type 
II outburst are not related to the orbital phase and imply an even 
higher increase in luminosity than Type I outburst s, reaching th e 
Eddington luminosity (for a review on BeXBs see lReigl(l201lh '). 

In the case of classical Supergiant X-ray Binaries (SGXBs), 
the compact object is embedded in the dense and powerful wind 
of a OB supergiant companion, swallowing everything that en¬ 
ters its gravitational domain. The mass loss rate of the donor 
is of the order of > 10^^ , and the compact object is 

usually found at a close distance of a ~ 1.5 - 2/?*. In such 
a close orbit, the captured matter is able to fuel a persistent 
X-ray emission of ~ 10^^“^® erg/s. Flares and off-states are 
often observed in SGXBs, indicating an abrupt transition in 
the accretion rate. They might be produced either by sudden 
variati ons of density in the medium transited by the compac t 
object (iMartlnez-Nilnez et al.ll2014t iKrevkenbohm et ani2008h . 
or either by instabilities above the magnetosphere of the neu- 
tron star, as propose d in the quasi-spherical accretion theory by 
IShakura et al.l (l2012h . 

The medium transited by the compact object through the ex¬ 
tended atmosphere of an OB supergiant star is not smooth be¬ 
cause of, at least, two phenomena. First, density inhomogeneities 


(clumps) are present as a n intrinsic feature of the radiatively 
driven winds of hot stars (iFucv & Whitel[T9^ lOskinova et al.l 
l2012h . Second, hydrodynamical simulations show that the X- 
ray radiation and the gravity field of the compact object disturb 
the wind of the donor, inducing the forma tion of denser struc¬ 
tures like filaments, bow shocks and wakes (iBlondin et alJll99n 

[mb- 

In the last decade and a half, new discoveries have led to the 
addition of new groups to the previous picture of HMXBs, stress¬ 
ing the value of grasping the different features of the sources 
such as geometry, compact object properties, optical star pecu¬ 
liarities and wind dumpiness. The new groups are Supergiant 
Fast X-ray Transient systems (SFXTs), y Cassiopeae analogs 
and y-ray Binaries. 

SFXTs are systems with a supergiant optical star, as in 
SGXBs, but defined by an extremely transient behaviour. During 
quiescence they exhibit low luminosity (~ 10^’^ erg/s), but they 
spend most of their time in an intermediate level of emission 
(~ 10^^^^"^ erg/s). They display short ou tbursts (~few hou rs), 
reaching luminosities up to 10^®^^^ erg/s (ISidoli et alj|2009h . It 
is likely that the dumpiness of the wind plays a main role in the 
variability of these sources. Other mechanisms involving cen¬ 
trifugal and magnetic barriers could enhance the observed lumi¬ 
nosity swings, relaxing the need ed variation ampli tudes in the 
physical conditions of the wind (iBozzo et alJl2008ll . Neverthe¬ 
less, other authors explain th e variability appea li ng to the quasi- 
spher ical accretion model (iDrave et al.l 120131 : iPaizis & Sidolil 

l2^ . 

y Cassiopeae analogs are characterized by the thermal nature 
of the X-ray emission, with plasma temperatures of ~ 10® K 
(~ 10 keV), an X-ray luminosity of 10®““®® erg/s, and high flux 
variability on various time scales. Howe ver, they do not display 
giant outbursts as observed in BeXBs (iFopes de Oliveira et al.l 
l2010h . Presently, it is not clear that the X-ray emission is emitted 
by accretion processes (onto a neutron star or a white dwarf), or 
alternatively generated from the interaction between the surface 
of the star, the circumstellar disk and its magnetic field. 

High Mass y-ray Binary systems (HMGBs) are HMXBs 
where the emission peaks above 1 MeV. Nowadays, it is thought 
that the emission is caused by accelerated particles in the shock 
that is produced when the pulsar wind collides the massive star 
wind. Therefore, they are powered by the rotational energy of 
the neutron star, in opposition to the rest of HMXBs, which are 
accretion fed. There are currently five confirmed HMGBs, all of 
the m with a main sequence optical star (for a review on HMGBs 
see lDubusl(l2()Tl 'l. 

Finally, there are sources which, following a number of rea¬ 
sons, cannot be classified in any of the already mentioned classes 
of HMXBs. Particularly, among the set of sources studied in this 
paper, they are 4U 2206H-54, Centaurus X-3 and Cygnu s X-1. 
The optical star in 4U 2206H-54 is a 09.5V dBlav et al.l 1200^ . 
neither a supergiant nor a Be star. The system may be part of 
a new group of w ind-fed HMXBs with a main sequence donor 
(iRibo et al.ll2006ll . Centaurus X-3 and Cygnus X-1 are the only 
systems here collect ed where accretion i s persistently drive n 
by an accretion disk (iTiemkes et alj[l9^ IShaniro et alJll97^ . 
what is reflected in the spectra of both sources. 

In this paper, we study the FeKa line for the whole sample 
of HMXBs available with XMM-Newton until August, 2013. In 
Section 12 we present the set of observations, the reduction pro¬ 
cess and the more important details concerning the spectral fits. 
In Section [2 we show the obtained results: an spectral atlas in¬ 
cluding every fit and different plots relating fit parameters. In 
Section 0] we interpret the obtained results and we summarize 
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the most important conclusions in Section |5] In Appendix A we 
present a set of tables describing the obtained parameters from 
the spectral fits. In Appendix B we show the spectral atlas, con¬ 
taining the plot of every spectrum that we have analysed in this 
survey. We show the observations and the models, together with 
the ratio between them. 

2. Observations and data treatment 

The XMM-Newton observatory (iLumb et al.ll2012h is fitted with 
three X-ray telescopes of 1500 crn^ and a coalignated optical 
telescope. Spectroscopy and photometry is done by the 6 instru¬ 
ments on board: three X-ray imaging cameras EPIC (European 
Photon Imaging Camera), two grating X ray spectrometers RGS 
(Reflection Grating Spectrometer) and an optical monitor (OM). 
EPIC cameras (0.1-15 keV) are the only instruments at XMM- 
Newton covering the energy range of the Ee complex. Among 
EPIC, one camera uses PN CCDs, and the other two use MOS 
CCDs. EPIC PN cameras (EPN) surpass by a factor ~3 the ef¬ 
fective area of the MOS cameras at 6-7 keV, making EPN more 
suitable for our purposes. Compared to other missions, the High 
Energy Transmission Grating Spectrometer (HETGS) on board 
of Chandra provides better energy resolution in the energy range 
of the Ee complex, but the effective area available with EPN is 
significantly higher. EPN provides the adequate conditions to 
perform the study here presented, on account of the moderate 
(but sufficient) spectral resolution (AEjE ~ 40), and great effec¬ 
tive area (~ 1000 cirP'), enabling us to analyze a large amount of 
sources in an homogeneous and consistent way. 

Since HMXBs are usually variable, we often observe in the 
same observation a dramatic change in luminosity, remarkably 
affecting the spectral parameters. In these cases, an averaged 
spectrum does not reproduce the actual emission of the source, 
and it is advisable to split the observation in more than one time 
interval. We have considered five different stated] of the systems 
in order to define the time intervals: dips, quiescence, flares, 
eclipse ingress/egress, and eclipse. We have used the follow¬ 
ing criteria. When luminosity drops a factor >2 in the timescale 
of <1 hour, we have tagged the time interval as a dip. Analo¬ 
gously, when luminosity rises >2 in the timescale of <1 hour, we 
have labelled the time interval as a flare. Eor observations cover¬ 
ing eclipsing phases, we have defined time intervals for eclipse 
ingress/egress and eclipse. The rest of time intervals are tagged 
as quiescent states. 

In Eigure [T] we see the light curve of an observation of 
4U 1538-522, as an example of how we have split the time in¬ 
tervals in the observations. The source was observed during the 
ingress in an X-ray eclipse, which is clearly noticeable in the 
light curve. We have separated the observation in two time inter¬ 
vals, one covering the ingress in eclipse, and another one cover¬ 
ing the eclipse. 

Summarizing the sample of observations, we have collected 
data from 46 HMXBs. 21 of them exhibit EeKa emission. We 
note that some sources have more than one available observa¬ 
tion. Taking everything into account (46 sources, temporal split¬ 
ting depending on the state of the source, and more than one 
observation per source in some cases), we end up with a total 
number of 108 spectra that we have a nalysed. _ 

We have followed the catalogue of iLiu et d1 (l2006l) . in addi¬ 
tion to later discoveries or confirmations, to identify the currently 
known HMXBs, and used every available XMM-Newton public 

* The states considered in this work and the also called states in black 
hole binary systems must not be mislead. 


Light curve 4U 1538-522 



Eig. 1: Light curve of the observation of 4U 1538-522 (Ob- 
sID:0152780201). We have split the observation in two parts, one for 
the ingress in the eclipse, and another one for the eclipse. 

observatioifl The sources not included in the Liu catalogue, but 
here considered, are: HD 119682, SS 397, IGR J16328-4726, 
HD 45314, HD 157832, Swift J045106.8-694803, IGR J16207- 
5129 and XTE J1743-363. 

2.1. Data reduction 

We have reduced the data using Science Analysis System (SAS), 
version 12.0.1. Since the sample of observations contain an het¬ 
erogenous group of HMXBs, we found different observation 
modes (timing and imaging) to account for the different prop¬ 
erties of the sources. In the brightest systems, the observations 
were usually performed using the timing mode, while the faintest 
sources were observed using imaging modes. 

Timing modes permit to process the arrival of photons at 
high rate, since only one CCD operates and the information is 
collapsed into one dimension, allowing a fast read out. The time 
resolu tion is as high as 30 fxs (7 fis in burst mode. iKrsch et all 
(l2006h L Even though the high timing resolution reached with 
these observation modes, pile-up is still present in several cases, 
specially when the count rate is > 800 counts s“'. We have 
checked in every observation if pile-up is affecting the data, us¬ 
ing the SAS task epatplot, and we have excised the core of the 
source’s point spread function in the pertinent cases. The size of 
the excised region has been chosen wide enough to remove the 
unwanted pile-u p effects (see examples of the use of epatplot in 
iNg etal.l(l2010h L 

Background subtraction process is also dependent on the 
brightness of the source. In the EPN timing mode, the PSE of 
the sources displaying > 200 counts s“' will span across the 
whole CCD. Therefore, any area selected as a background re¬ 
gion will be contaminated by source photons. Since this effect 
is strongly energy dependent, for the brightest sources we have 
chosen a method of background subtract ion similar to the one 
perfor med in the analysis of Vela X-1 by iMartlnez-Nunez et al.l 
(l2014h . where a blank sky spectrum taken in timing mode is used 
as the real background for energies below 2.5 keV, while the rest 
of the spectrum corresponds to the outermost pixels of the CCD. 

^ http://xmm.esac.esa.int/xsa/ 
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Reduced-^^ 


Fig. 2: Number of accepted models depending on the Reduced-;^'^ 
value. 


Meanwhile, for common observations, we have used source-free 
regions to extract a background spectrum and subtract it from 
the former source plus background energy distribution. 

Ancillary response files were generated using the SAS task 
arfgen. For observations taken in timing mode affected by pile- 
up, we have followed the recommendations of the XMM-Newton 
SAS User Guide in order to generate the appropriate ancillary 
response files. Response matrices were created using the SAS 
task nnfgen. 

2.2. Spectral fitting 

For the spectral analysis we have used XSPEC, version 12.8.cB 
We have rebinned the spectra to have a minimum of 20 counts 
per bin, and a bin size of at least 1/3 of the FWHM of the intrinsic 
energy resolution, in order to be allo wed to apply statistics in 
the fitting of a set of poissonian data (ICas i^ll979h . 

In Table [T] we present the sample of models employed for 
the continuum in the fits. Every model is a combination of ad¬ 
ditive and multiplicative models. An additive model stands for 
a source of X-rays (e.g. bremsstrahlung radiation), and a multi¬ 
plicative model represents a energy-dependent change of an ad¬ 
ditive model (e.g. photoelectric absorption). 

The models presented in Table [1] have been tested in every 
observation, and accepted depending on the Reduced-;^^^ 
with n the number of bins and m the number of fitted parame¬ 
ters). Every observation has particular characteristics, and there¬ 
fore the decision of which Reduced-;^^^ value is acceptable has 
been taken one by one. In Eigure|2]we can see that most of the 
fits result in a Reduced-;^^^ ^ 1, as expected for a suitable fit. 
The highest value of Reduced-;^^” for an accepted model has been 
1.82. The parameters arising from the fits are listed in Tables lATSl 
and lA.3l 

We can classify the additive components of the models as 
thermal or non-thermal. A component is called thermal when 
radiation is produced as a consequence of the thermal motion 
of the plasma particles (e.g. blackbody radiation). Otherwise, 
the emitted radiation is non-thermal (e.g. non-thermal Inverse- 
Compton emission). If all the additive components of a model 

^ http://heasarc.nasa.gov/xanadu/xspec/ 


are thermal, we classify the model as thermal (analogously for 
non-thermal). We have also used hybrid models, combining ther¬ 
mal and non-thermal components. The thermal components used 
in this work are the following: 

- bbody: blackbody emission. 

- diskbb: model of an accretion disk emission made of multiple 
blackbody components. 

- bremss: thermal bremsstrahlung emission (electrons dis¬ 
tributed according to the Maxwell-Boltzmann distribution). 

- mekal: emission from optically thin hot gas, includ ing spec¬ 
tral lines from several elements (iMewe et al.lll985h . 

- cemekl: built from the mekal model, incorporating multi¬ 
temperature emission. 

On the other hand, the only non-thermal component used in 
this work is: 


- powerlaw: phenomenological model consisting of a simple 
inverse power law profile (oc E~^). This profile is a foot¬ 
print of Inverse-Compton scattering by hot electrons (non- 
thermally distributed) of a seed radiation field. 


Eor the photoelectri c absorption , we h ave used tbne\^ the 
improved version of thetWlrns,e^^ (l2000h model tbabs, setting 
the cross sections t o the Verner et al. jll996h ones and the abun¬ 
dances according to lWilms et alJ ( 20001) . The most important pa¬ 
rameter of this model is the total equivalent hydrogen column 
Nh, which is the integrated amount of hydrogen atoms in the 
line of sight from the observer to the source, per cm^. We have 
also added the model cabs to account for the Compton scatter¬ 
ing, which is not comprised in the tbnew model and is especially 
significative for A// > 10^"^ cmT^. 

The emission lines are fitted using Gaussian profiles. We 
have categorized as EeKcr any emission line that fulfils the fol¬ 
lowing conditions: 


1) The centroid energy of the Gaussian component lies in the in¬ 
terval [6.3,6.65] keV. The interval includes the expected en¬ 
ergy of EeKa emission from Ee II ( ~ 6.395 keV) to Ee XXIII 
(~ 6.63 keV) dKallman et alJl2()()4l) . This condition excludes 
the detection of any hypothetical fluorescent emission from 
Ee XXIV-XXV at ~ 6.67-6.7 keV, thereby excluding any con¬ 
fusion between EeKa and the recombination line Ee XXV at 
similar wavelength. The fluorescence yields of Ee XXIV-XXV 
are low compared to lower ionization states. 

2) The statistical significance (cr^j^,,) of the Gaussian com¬ 
ponent is greater than 2cr. We have calculated cTsign from 

xl, -xlp assuming ^2 where 

from a fiit using certain model with the Gaussian component 
included, and xl^ arising from a fit using the same model 
without the Gaussian component. 

In some cases, EeKa line is clearly noticeable, but EeKyS is 
not prominent enough to permit error estimation of its parame¬ 
ters. In these cases, we have co nstrained the centroid energy and 
the norm of EeK/I according to lKallman et alJ (l2004l) : 


- Energy(EeKj6) = Energy(EeKa) + 0.652 keV 

- Norm(EeKj6) = Norm(EeKa) x 0.13 photons/cm^/s 

2 http://pulsar.sternwarte.uni-erlangen.de/wilms/research/tbabs/ 

^ This assumption is not strictly true, since andxf., are not indepen¬ 
dent. However, it provides an estimation on the impact of the Gaussian 
component in the model. 
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The estimated parameters, like the EW, are very sensitive to 
the fit of the continuum. Therefore, although the Fe complex 
appears in the ~6-7 keV energies, we broadened the spectral 
scope to an energy range of 1-10 keV to perform the analysis. 
It also allows us to consider possible calibration inaccuracies 
in the charge transfer inefficiency (CTI) and the X-ray loading 
(XRL), a n issue reported in previous a naly sis of EPN observ a- 
tions (see lMartmez-Nunez et al.l (1201 4l) an d lFlirst et al.1 J^Onl) '). 
In the few cases of possible CTI or XRL, we applied an artificial 
gain E„ew = + offset (see Table|A.2t. 

The estimation of the parameter confidence regions (at 90% 
level) have been calculated with a Markov Chain Monte Carlo 
(MCMC) technique, implemented in XSPEC, where N genera¬ 
tions of the set of free parameters are used to determine the best- 
fit values and the confidence regions. We have set A = 1.5 x 10^ 
in our calculations. These chains are also valid to estimate fluxes 
and equivalent widths. 

3. Results 

3.1. Spectral Atlas 

In Appendix lEl we present the full sample of analysed spectra. 
The figures in AppendixiBishow the set of analysed observations 
(cross points), the model employed (solid line), the components 
of the model (dotted line) and the ratio between observation and 
model (lower box in each spectrum plot). 

We show a list of the sources in Table lA.ll giving the class 
where we have grouped them and the reference for such a clas¬ 
sification. We can see that the different classes of HMXBs be¬ 
have qualitatively different in the region of the Fe complex (~6- 
7 keV), reflecting the distinct accretion regimes that characterize 
them. We have observed three patterns in the Fe complex, that 
we define as Type I, II and III (see Figure^. We define Type I, 
when fluorescence lines FeKa and FeKyS are observed, but not 
recombination lines Fe XXV and Fe XXVI. We define Type II, 
when fluorescence lines are detected, together with recombina¬ 
tion lines Fe XXV and Fe XXVI. Finally, we define Type III, when 
Fe lines are not detected. 


The general features observed in this work, for the different 
groups of HMXB, are summarized in Table|2] and explained be¬ 
low in more details: 

- BeXBs. We have collected data from 10 sources. All the ob¬ 
servations were performed in quiescence. We have detected 
FeKa emission in only one BeXB (SAX J2103.5+4545). The 
upper limit of the FeKa EW in the rest of BeXBs is in gen¬ 
eral higher than the observed value in SAX J2103.5+4545, 
implying that the lack of detections might be due to a poor 
signal-to-noise. The spectra can be modeled by thermal or 
a combination of thermal and non-thermal components, ex¬ 
cept for Swift J045106.8-694803 (fitted using an absorbed 
power law). Seven sources accept a thermal model, and 6 a 
combination model (4 of them accept both). 

- SGXBs. We have gathered data from 12 sources. Ten of them 
show detectable Fe fluorescence emission. The only excep¬ 
tions are IGR J16465-4507 and SAX J1802.7-2017,the most 
distant SGXBs included in this work, at 12.5 and 12.4 kpc re¬ 
spectively. The EW upper limits in these two sources is high, 
implying that their faintness is very likely the reason why we 


Group 

# Sources 

Fe complex 

Models 

Nh 

BeXB 

10 

Type III 

T, TN 

Low 

SGXB 

12 

Type I 

N 

High 

SFXT 

10 

Type III 
(in quiescence) 

T, N, TN 

High 

y Cass like 

8 

Type II 

T 

Low 

HMGB 

2 

Type III 

T, N, TN 

Low 


Table 2: Description of the features observed in this work, for the 
different groups of HMXBs. We have analyzed data of 46 sources. 
However, those classified as peculiars (3 sources) or non classified 
(AX J1749.1-2733) are not included in this table. We define Nh of a 
group as high, when the typical value observed is well over the estima- 
tions of the interstellar in the line of sight of the sources following 
IWillingale et all (l201.3ll . That is, we say that the Nh of a group is high 
when the absorption is typically intrinsic of the systems. 

do not detect FeKa. The 12 SGXBs can be well fitted using 
non-thermal models, although thermal components are also 
plausible in some sources. In general, SGXBs are character¬ 
ized by high absorption and the presence of Fe fluorescence 
emission lines. 

- SFXTs. We have collected data from 10 sources. Three of 
them show FeKa: AX J1841.0-0536, IGR J11215-5952 and 
IGR J16479-4514. The EW upper limit in the rest of sources 
is high. Therefore, FeKa would be probably detectable with 
a better signal-to-noise. The models employed for fitting the 
SFXT systems are very heterogeneous, with no preference 
of thermal, non-thermal, neither a combination of both kinds 
of models. 

- y Cassiopeae analogs. We have gathered observations from 
8 sources. Five of them exhibit FeKa. The EW upper limit 
in the other 3 sources is very high. Again, it implies a 
very likely presence of fluorescence in the case of better 
signal-to-noise. In addition, recombination lines of Fe XXV 
and Fe XXVI, are always present in the set of y Cassio¬ 
peae analogs. These lines are included in the XSPEC model 
mekal. For most of the observations we have achieved a 
good fit using a combination of mekal components. In a few 
cases we have used other components: diskbb and power- 
law, but mekal is by far the most employed one in y Cassio- 
peae like systems, in agreement with previous X-ray analysis 
(iLopes de Oliveira et al.ll2010ll200^ . 

- HMGBs. We have collected data from two HMGBs: 
LS 1+61 303 and LS 5039. None of them show Fe features. 
However, the signal-to-noise in these observations is poor 
and the upper limits of the FeKa EW are high enough to 
do not rule out the presence of the line. We have used both 
thermal and non-thermal components in the fits. 

- Peculiars. Set of sources that do not accommodate in any 
of the aforementioned classes of HMXBs, as explained in 
the introduction. We have collected data of 3 such systems: 
4U 2206+54, Centaurus X-3 and Cygnus X-1: 

- 4U 2206+54 does not show any detectable Fe emission 
line, and the upper limit in the EW of FeKa is low (it is 
comparable to the upper limits in the BeXBs). It can be 
fitted by means of an hybrid model (thermal plus non- 
thermal components). 
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Model Continuum models 


Non-thermal 

N2 

Ni 

powerlaw x tbnew x cabs 

powerlawi x tbnew\ x cabs\ + powerlaw 2 x tbnew 2 x cabs 2 
powerlawi X tbnew\ x cabs\ + powerlaw 2 x tbnew 2 x cabs 2 + powerlaw^, x tbnewj, x cabsj, 


Ti 

mekal X tbnew x cabs 


Ti 

(mekal -I- mekal) X tbnew X cabs 


T3 

(mekal + mekal + mekal) X tbnew X cabs 


Ta 

(cemekl) X tbnew X cabs 


Ts 

bbody X tbnew x cabs 


Te 

(bbody + bbody) x tbnew x cabs 

Thermal 

Ti 

(bbodyi -H bbody 2 ) x tbnew x cabs 


T& 

bbodyx x tbnew\ x cabs\ -i- bbody 2 x tbnew 2 x cabs 2 


Tg 

diskbb X tbnew x cabs 


Tio 

(diskbb + bbody) x tbnew x cabs 


Tn 

bremss x tbnew x cabs 


Tn 

bbody X tbnew x cabs + bremss x tbnew x cabs 


Tii 

(bremss + bbody) x tbnew x cabs 


TNi 

(powerlaw + bbody) x tbnew x cabs 


TN2 

powerlaw X tbnewi x cabsi + bbody x tbnew 2 x cabs 2 

Both 

TNs 

(powerlaw + diskbb) x tbnew x cabs 


TNa 

powerlaw X tbnewi X cabs\ + diskbb x tbnew 2 x cabs 2 


TNs 

(powerlaw -I- mekal) X tbnew X cabs 


Table 1; List of models used to fit the continuum, described in XSPEC notation. The basic components are powerlaw, bbody, diskbb, breniss, 
mekal and cemekl, together with tbnew to account for the absorption and cabs for the non-relativistic Compton scattering. The employed models 
are a combination of these components, in addition to Gaussian profiles modelling emission lines. We divide the models in three types: for 

non-thermal, T# for thermal and TN# for models containing both thermal and non-thermal components. 


- Centaurus X-3 presents a rich emission lines spectrum. 
Concretely, in the Fe complex we are able to identify 
FeKa, FeK/l, Fe XXV and Fe XXVI. We have used either 
an hybrid model either a non-thermal model. 

- Cygnus X-1 exhibits a broad Fe feature, sometimes com¬ 
bined with a faint and narrow, but statistically significant, 
FeKa line. We have mostly used non-thermal models, oc¬ 
casionally combined with a bbody or diskbb component. 

- AX J1749.1-2733. In thi s system, the optica l member has 
been classified as a Bl-2 dKarasev et alJl201(ih . but the lumi¬ 
nosity class remains unknown, preventing us to incorporate 
the source in a defined group. Although it does not exhibit 
detectable Fe emission, the high absorption clearly points to 
a supergiant companion. In addition, the EW upper limit of 
FeKa is compatible with the values observed in SGXBs and 
SFXTs. It can be well fitted using an absorbed powerlaw or 
a blackbody. 

In summary, it is very likely that FeKa is an ubiquitous fea¬ 
ture in HMXBs, and its detection highly depends on the qual¬ 
ity of the observations. In this regard, the EW of the line 
is very affected by the level of intrinsic absorption present 
on the sources (see also Section [3.4.3l) . SGXBs and SEXTs, 
which show higher absorption, tend to exhibit a more promi¬ 
nent Ee fluorescence. 


3.2. FeKa width 

In Table IA.3I we show the parameters of every detected EeKa, 
including the width of the line (cr;,„g). We have made a dis¬ 
tinction between narrow and broad EeKa. We have defined 
narrow EeKa when crune < 0.15 keV, and broad EeKa when 
crime > 0.15 keV. This separation is both physically and observa- 
tionally motivated: 


The origin of broad Ee features in X-ray Binaries is still an 
opened question, but the most likely alter natives are re l ated t o 
the presence of an accretion disk (see e.g. iHanke et~^ (l2009l) : 
iNg alJ (l2010h : iDuro et ^ (1201 ih l. However, narrow features 
are not compatible with material rotating at high velocities or rel- 
ativistically broadened. Given that broad and narrow EeKa have 
a clearly different origin, they must be analysed in a separated 
way. Then, it raises the question of how to define the separation 
mark between them. In Eigure|4]we can see that in our sample 
the number of detected sources decreases when increasing cr;,„e. 
Moreover, most of the detections are grouped at cr;;„g <0.15 keV. 
Hence, cr;,„g <0.15 keV seems a natural threshold for the def¬ 
inition of narrow lines in the sample. In addition, we must pay 
attention to the plausible contamination of EeKa with Ee XXV, 
which is located at ~ 6.7 keV. The chosen criterion separates the 
sources where it is very unlikely that EeKa is contaminated by 
Ee XXV (narrow lines), from the sources which probably suffer 
from this issue (broad lines). A more detailed analysis of broad 
Ee features in HMXBs is out of the scope of this paper and its 
discussion will be fully developed in a forthcoming work. 

Hereafter, when FeKa is mentioned, we refer to the narrow 
feature. Erom the total number of 108 analysed observations we 
have detected (narrow) EeKa in 60 of them. 

3.3. Centroid energy 

In Eigure|5]we can see an histogram presenting the centroid en¬ 
ergy of EeKa. A Gaussian fit of the data reveals a mean value 
for the centroid energy of 6.42 keV. There are not significant 
differences in the averaged values obtained neither for different 
classes of HMXBs, neither for different states. The standard de¬ 
viation is 0.02 keV, comparable to the the error that we typically 
obtain in the estimation of the centroid energy in the fits (see 
Table IAJTi. Taking into account the standard deviation and the 
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Fig. 3: Patterns found in the Fe complex of HMXBs: Type I (left panel). Type II (central panel), and Type III (right panel). 
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Fig. 4; Histogram of the FeKa width. The bulk of the detections are 
grouped showing (T;„„ < 0.15 keV. We define them as narrow FeKa. 
The rest are defined as broad Fe features. Even though we have detected 
60 narrow FeKa, this plot only includes 38. The reason is that 22 of 
them are very narrow (or the signal-to-noise very low) and they have 
been treated in the fits as delta functions (in Table lA.3l we present their 
width as = 0). 


Fig. 5: Centroid energy of FeKa, with a Gaussian fit overplotted (blue 
profile). The mean value is 6.42 ± 0.02 keV, compatible with Fe I-XVIII. 
Even though we have 60 detections of FeKa, in this plot we only see 
55. Five cases fulfil the requirements of detection, but the low signal-to- 
noise do not permit to find an accurate centroid energy. They have not 
been included in this plot. In these five cases we have set the centroid 
energy of FeKa to ~6.4 keV. 


uncertainties in the CTI corrections in EPhQ, the centroid en¬ 
ergy of FeKa constrains the ionizat ion state of Fe to less ionized 
than Fe XVIII (iKallman et al.||2004 ), in agreement with previous 


studies in HM XBs (iTorreionet al.ll2010bt Gottwalde^l. 19^ 


I 


lNagaselfT989h . In this regard, the study of ITorreion et al.l (12010 
using HETGS (more accurate in wavelength than EPN), gives a 
narrower constrain in the ionization state (Ee l-x). Our present 
work support this result adding more sources to the sample. 


In the right side of Eigure|5]we can see seven EeKa detec¬ 
tions emerging out of the Gaussian profile. Eour of them are un¬ 
likely to be described by such a Gaussian profile, since they lie 
more than three times the standard deviation away from the mean 
energy. All four belong to Cygnus X-1. 


® Please find more information about long-term CTI correction in the 
release note EPlC-pn Long-Term CTI, by M.J.S Smith et al. (2014), at 
http://xmm2.esac.esa.int/docs/documents/CAL-SRN-0309-l-0.ps.gz 
and EPIC status of calibration and data analysis by Guainazzi (2008), 
at http://xmm2.esac.esa.int/extemal/xmm_sw_cal/calib/CAL-TN- 
0018.pdf. 


3.4. Correlated parameters 

One of the goals of this work is to study plausible correlations 
involving the parameters of EeKa (position, width, intensity and 
EW, and other parameters in the fits, such as the absorbing col¬ 
umn and the intensity of the continuum. Note that, even when a 
good fit is reached, the confidence region of a parameter might 
be occasionally difficult to find due to the dependence of the pa¬ 
rameter with other parameters of the model. We specify in each 
of the following subsections the number of cases where a suc¬ 
cessful estimation of the 90% confidence region has been done. 

3.4.1. Continuum Flux vs FeKa Flux 

In Eigure|6]we have represented the unabsorbed flux of the con¬ 
tinuum between 1-10 keV cancelling EeKa emission (Fi_iojtcy), 
against the flux of EeKa (FfcKa)- We have successfully found a 
90% confidence region of the flux of EeKa in 56 cases. 

In a logarithmic scale, we identify two different patterns of 
correlation. Eirst, for a subset including all the eclipse observa¬ 
tions and IGR J16318-4848, we And a correlation with Pearson 
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FeKa Flux (erg/s/cm^) Continuum Luminosityi_^Qj^gv (erg/s) 


Fig. 6: Fi-iohev versus FfeKa- Blue dashed line marks the correlation 
observed for IGR J16318-4848 jointly with eclipse observations, and 
the red solid line follow the bulk of the observations. The color map 
indicates the (T„j„ of the line (defined in Section [2^ . 


Coefficient (PC) of 0.98 (blue dashed line). Second, for the rest 
of observations, we find a correlation with PC=0.89 (red solid 
line). 

A linear fit of the parameters (in logarithmic scale) in the out- 
of-eclipse observations (red line in Figure|6l), gives the following 
dependence; 


Fi-iofev(erg/s/cmO = X 10 


-v2.18±0.87 


( 1 ) 


The errors show the standard deviation of the parameters in the 
linear fit. 

The observed divergence among eclipse (plus IGR J16318- 
4848) and out-of-eclipse observations suggests that the compan¬ 
ion star blocks in a different proportion the continuum and the 
FeKcr emission. Therefore, an important contribution of the fluo¬ 
rescence emission is produced in an extended region of R > R*. 
This is consistent w ith previous analysis of ecli pse observa- 
tions of HMXBs (e.g. Rodgs-Roca_etalj( 201 Ih and lAudlev et alJ 
(120061) '). In particular, Audlev et alJ ( 20()6h estimates that 20% of 
FeKcr in OAO 1657-415 is emitted from 19 light-seconds off the 
X-ray source. 

We have also the luminosity of the continuum, in order to 
face it with the EW of FeKa. For the flux-to-luminosity conver¬ 
sion we have used the estimations of the distance shown in Ta- 
ble lA.ll We have excluded eclipse and IGR J16318-4848 from 
this analysis, given that the EW is strongly affected by the high 
obscuration of the continuum that they suffer from. In Eigure Q 
we plot the EW of EeKa against the unabsorbed luminosity 
of the continuum between 1-10 keV cancelling EeKa emission 
(Li-iOitcv). We observe two different groups of sources: 1) y Cas- 
siopeae analogs lie at low luminosities (Li-iofev < 10^^ erg/s); 
2) the rest of sources which exhibit EeKa. y Cassiopeae analogs 
don’t show any evident correlation (they are very few points), 
while the rest present a moderate inverse correlation (PC=-0.25, 
and PC=-0.39 using only the sources with an available estima¬ 
tion of distance with error, marked as filled diamonds in Eig¬ 
ure |7](. A linear fit using the filled diamonds in Eigure |7] gives: 




(2) 


Eig. 7: EW of FeKo- against Li-iotcy. y Cassiopeae analogs (circles) 
lie at < 10'^^ erg/s. Open symbols indicate that either the dis¬ 

tance either the error in the estimation of the distance is unknown. The 
solid line corresponds to a linear fit in logarithmic scale of the filled dia¬ 
monds, that is, the sources with available distance with error estimation 
and > 10^^ erg/s. 


iBaldwinI (Il977h observed an inverse correlation in the EW 
of CW and the LTV luminosity in AGNs. Analogously, the de¬ 
crease of the EW of EeKa when increasing the X-ray luminos¬ 
ity is called X-rays Baldwin effect. The dependence that we ob- 
serve is comp a tible w ithin the error with the one observed by 
iTorreion et al.l (1201 Obh in X-ray Binaries using a narrower en¬ 
ergy range: EW « L-O-29 

3.4.2. FeKa Width vs Centroid Energy 

In Eigure [8] we present the centroid energy of this feature ver¬ 
sus its width {(Jiine)- We have successfully found a 90% confi¬ 
dence region of cr;,„e in 20 cases. We can see a moderate cor¬ 
relation (black line in Eigure PC=0.55), indicating a pos¬ 
sible blending of lines. Two observations (uppermost side of 
Eig. ISJ do not follow the correlation. They correspond to ob¬ 
servations of 4U 1700-37 (Obs.ID 0083280201) and EXO 1722- 
363 (Obs.ID 0405640201) where the Ee complex is hardly re¬ 
solved, and therefore it is very likely that a contribution of 
Ee XXV and Ee XXVI in the model of EeKa is increasing the 
measured width of the EeKa line. 

Coloured squares correspond to the expected width from 
the contribution of three different broadening phenomena: line 
blending, Doppler shifts and Compton broadening. A discussion 
on the different broadening contributions is given in SectionH] 

3.4.3. Curve of growth 

In Eigurel^we show, for out of eclipse observations, the Nh ver¬ 
sus the EW of EeKa (what is generally known as the curve of 
growth). We have successfully found a 90% confidence region 
of both Nh and EW in 46 cases. We want to take into account 
observations were Nh reflects the intrinsic absorption of the sys¬ 
tem. Hence, we have set Nh > 2 as a condition to safely ex¬ 
ceed the typical Nh of the interstellar medium for the sources 
here studied (checked using the online application following 
IWillingale et al.l (1201 3h l. The use of this criterion excludes the 
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FeKa Centroid Energy (keV) 


Fig. 8: Width of FeKa icrii„g) versus the centroid energy (black 
squares). Black solid line traces a linear fit. We have marked in 
colour the expected width from the contribution of three differ¬ 
ent broadening phenomena: line blending, Compton broadening and 
Doppler shifts, considering velocities of V(km/s) = 1000 (red) and 
V(km/s) = 2000 (green). Every single black square has associated a sin¬ 
gle red square and a single green square corresponding to the expected 
values of ct/,-,,,, for that observation (see Section|4j. 

BeXB SAX J2103.5H-4545, the y Cassiopeae analogs: y Cassio- 
peae and HD 110432; and the SFXT IGR J11215-5952. More¬ 
over, eclipse observations show higher EW and they are not com¬ 
parable to out-of-eclipse observations. Therefore, eclipse obser¬ 
vations have not been plotted in Figure 0 As a consequence of 
the chosen criteria, we end up with a set of 36 observations, 
where all the donors are supergiants. 

Nh and the EW of EeKa are expect ed to correlate in 
HMXBs, as shown bv lTorreion et al.l (1201 Obi) , since the spectral 
lines are usually stronger when the optical depth increases. Our 
sample confirms these expectations, showing a notable correla¬ 
tion (PC=0.85). 

We have determined the theoretical curve of growth using 
numerical simulations. In this simulations there is an input of X- 
ray radiation with a power law profile, that is tra nsmitted throug h 
a cloud of spherically distributed neutral matter (lEikmannl2012l) . 

We have taken into account the power law index (F) in the 
simulations, since steeper profiles entail less photons available 
above the Fe K-shell threshold energy, decreasing the EW. That 
is, for the same Nh, the higher is F, the lower is the EW of EeKa. 
In Eigurel^the turquoise band traces the results from the simula¬ 
tions with F 6 [0.5,2], which is the typical range where we find 
F in our fits. 

3.5. Nh: SGXBs and SFXTs 

In Figure [TO] we have plotted histograms for the Nh values ob¬ 
served in SGXBs (red) and SFXTs (blue). In the cases where we 
have more than one observation for the same source, we have 
averaged the values in order to obtain one Nh representative for 
each system. The orbital phase critically affects the observed Nh, 
and therefore ingress/egress and eclipse phases have not been 
taken into account. 

We find that SGXBs are in general more absorbed than 
SFXTs. We have performed a permutation test in order to quan¬ 
tify if the observed disparity in the Nh is a random effect. We 



Nequ'v^(102z atoms cm 


Fig. 9: Curve of growth observed for FeKa, that is, EW against 
Nh. The turquoise band marks the expected correlation using nu¬ 
merical simulations. The sources are identified by different symbols 
when more than one observation is included: 4U 1700-37 (open cir¬ 
cle), 4U 1907-1-09 (open upward triangle), Cygnus X-1 (open down¬ 
ward triangle), EX01722-363 (open diamond), IGR J16318-4848 (open 
square) and IGR J16320-4751 (plus). Only one observation for Cen- 
taurus X-3, GX 301-2, Vela X-1 and XTE J0421-I-560 (all four a star 
symbol). 
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Fig. 10: Histograms showing a comparison of the Nh values observed 
in SGXBs (filled red) and SFXTs (empty blue). 


have 10 Ah values for SGXBs and 9 for SFXTs. We have merged 
them in a set of 19 elements, and considered every possible com¬ 
bination of 2 groups of 10 and 9 elements (92378 possibilities). 
We have compared the median of the two subsets, and calcu¬ 
lated the absolute difference. 99.7% of the cases have produced a 
lower absolute difference than the observed one. Using the mean 
instead of the median the percentage is also very high (98.8%). 
In conclusion, it is very likely that the discrepancies in the ob¬ 
served Nh values for SGXBs and SFXTs are produced by phys¬ 
ical reasons rather than they arise by chance. 
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Fig. 11: EW of FeKo- against Nh, in IGR J16320-4751. The turquoise 
band marks the numerical simulations results, with T e [0.5, 2], 


3.6. Individual sources analysis: IGR J16320-4751 and 
4U 1700-37 

3.6.1. IGR J16320-4751 


IGR J16320-4751 was detected by ASCA in 1994 and 1997 (cor- 
responding to A X J1631.9-4752), and by INTEGRAL in 2003 
(iTomsick et al.ll2()()^ . It is a HMXB composed b y an 081 op¬ 
tical star and a neutron star (iRahoui et al.l 1200^ . It shows a 
modu lation of 8.96 day s, that it is considered its orbital pe- 


riod ( Corbeteta. 20051) . and a pulsation period of ~ 1300s 


(iLutovinov et al.l 2005h . The ESA archives permits to collect 
eleven observations of IGR J16320-4751, enabling us to study 
in more detail the curve of growth, as well as tracking the ab¬ 
sorption variation along the orbital phase. 

In Fig. m we can see the curve of growth, as shown in 
Fig. 12 restricted to IGR J16320-4751. We clearly see the de¬ 
pendence between Nh and EW of FeKcr, as stated for the bulk 
of the sources in Section 13.4.31 and the general trend follow¬ 
ing the numerical simulations marked with the turquoise band. 
However, the agreement with the simulations is not completely 
fulfilled, given that the spectral fits of IGR J16320-4751 have 
brought out a power law index F ~ 0.5 (see Table Ia!^ . Since 
from a lower power law index we expect more EW of FeKcr, the 
points for IGR J16320-4751 are expected to be located in the up¬ 
per edge of the turquoise band, corresponding to the simulations 
with F = 0.5. We consider that the general trend is correct, but 
there are still some uncertainties in the fits and/or the theoretical 
hypothesis (spherical geometry and neutral matter). 

From 14th August to 17th September, in 2008, a campaign 
consisting of nine observations of IGR J16320-4751 was per¬ 
formed by XMM-Newton. We have used this set of data to plot 
the Nh modulation depending on the orbital phase (Fig.fTSli. We 
have set 0 = 0 at the Nh maximum. We have also calculated the 
theoretical absorption expected from a smooth wind in a non ec¬ 
centric orbit using a fl velocity law (ICastor et al.ll 19751) and the 
motion equation, taking into account variations in the orbital in¬ 
clination /, orbital separation a, star radius R*, mass loss rate M, 
parameter yS and the terminal velocity of the wind Voo. 

Indeed, M and Voo cannot be disentangled in this simple 
model, and the actual parameter used is M/voa. However, here¬ 
after we give values of M and Voo as if they were free vari¬ 


Fig. 12: Orbital modulation of Nh in the system IGR J16320-4751. 
Solid lines correspond to the expected absorption from a smooth wind 
and a non eccentric orbit, assuming an orbital separation a = 1.6 R*, 
R* = 20Rq, M = 10“^ Mglyr, v„ = lOOkmIs, p = 0.8; and different 
orbital inclinations i = 0, ;r/10, 7t/6 rad (red, green and blue). 


ables, since they are much more commonly used than M/uoo 
in the literature. This way, we constrain our parameters to the 
observed and predicted range of values in O supergiants: M = 
10“^ - -10~^ Mptly i and Voo = 500 - 3000 km/s (iK^udritzki & Pulsl 
l2000l:IVink et~^l2001h . 

For a null orbital inclination, we obviously obtain a flat Nh 
modulation (red line in Fig. fTSI) . that describes the observed Nh 
(except at (p = 0), assuming a - 1.6 R*, R* = 20Ro, M = 
lO^^Mo/yr, Uoo - 700 km/s and fi - 0.8. Gradually increasing 
the orbital inclination (/ = tt/IO, njb rad; green and blue lines 
respectively), we are able to describe a high Nh at 0 = 0, but 
loosing similarities around other orbital phases. Anyhow, given 
the simplicity of the model, the obtained parameters are certainly 
just indicatives. 


3.6.2. 4U 1700-37 

4U 1700-37 was detected for the first time by Uhuru in 1970 
dJones et aDll973h . The optic al star is HD 1539 19, a 06.5Iaf lo¬ 
cated at a distance of 1.9 kpc (lAnkav et al.l200lh . The orbital pe¬ 
riod is 3.41 days. Since X-ray pulsations have not been detected 
so far, the compact object can be either a neutron star or either 
a black hole. The database of ESA contains five observations of 
4U 1700-37, that we split in nine spectra in order to distinguish 
different states of the source. 

In FigurefOlwe can see the curve of growth, for 4U 1700-37. 
Although seven of the nine spectra show FeKa, we were able to 
constrain the boundaries of Nh in only five of the analysis. One 
of them corresponds to an eclipse observation (filled circle). It 
shows much more EW because the continuum flux is blocked by 
the optical star, whereas FeKo- comes from a more extended re¬ 
gion that is not completely hidden during eclipse. We do not see 
an obvious dependence between Nh and EW, although the points 
lie in a region close to the expected values (turquoise band). Any¬ 
way, a set of four observations (excluding the eclipse), is too 
small to perform an statistical analysis. 

From 17th to 20th of February 2001, 4U 1700-37 was ob¬ 
served by XMM-Newton four times, in a campaign covering dif- 
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Fig. 13: EW of FeKo- against Nh, in 4U 1700-37. The filled circle 
corresponds to an eclipse observation. The turquoise band traces the 
numerical calculations with T 6 [0.5,2]. 



Fig. 14: Orbital modulation of Nh in 4U 1700-37. Solid line correspond 
to the expected absorption from a smooth wind and a non ec centric or¬ 
bit, assuming the stellar values obtained in IClark et all ll2002h . an orbital 
inclination i ~ nil rad, and an orbital separation a = 1.41?*. 


ferent orbital phases. Therefore we can study the orbital modu¬ 
lation of Nh in the same way as we have done with IGR J16320- 
4751, but i ncluding more con strains coming from the non-LTE 
analysis of IClark et alJ (l2002h . where the following parameters 
are derived: /?* = 21.9Ro, M - 9.5 x 10“^Mo/yr, Voc - 
\150kml s and jS - 1.3. Taking into account that it is an eclips¬ 
ing binary, we assume i ~ njl. Therefore the only free parame¬ 
ter in our toy model is the orbital separation a. The best agree¬ 
ment is achieved when a - 1.4/?* (see Figure O. This orbital 
separati on is consistent (in abs olute units) to previous estima¬ 
tions of Conti^_Cowde^ (Il975h (/?* = 20 Rq, a - 1.35/?*) and 
iHeap & Corcoranf i 1992 ) (/?* = 18 + 3 /?©, a = 2.0 ± 0.4/?*). 


4. Discussion 

In Figure |5] we have shown the centroid energies of FeKa. The 
distribution of the histogram is roughly Gaussian, with a stan¬ 


dard deviation that reflects the uncertainties in the fits. However, 
four values are too high to be compatible with this distribution, 
and all of them belong to Cygnus X-1. It can be caused either 
by a not adequate fit or either by a physical reason. As stated be¬ 
fore, in Cygnus X-1 we detect a broad Fe feature, interpreted as 
a relativistically broadened fluorescence line. However, we have 
modelled the relativistically broadened feature with a Gaussian 
profile, that gives an acceptable fit, but might be inadequate, af¬ 
fecting the parameters of the narrow FeKa arising from the fits. 
Alternatively, as a plausible physical explanation, in Cygnus X-1 
the matter is accreted via an accretion disk, in opposition to the 
wind-fed accretion of most of the sources showing FeKa in this 
study. Therefore, the physical properties of the region emitting 
fluorescence might be different in Cygnus X-1 and the rest of 
systems. If this region is hotter in Cygnus X-1, the centroid en¬ 
ergy of FeKa would be shifted to higher energies as we observe. 

Regarding IGR J16318-4848, it is one of the most absorbed 
systems of the sample, and presents a special configuration of 
matter in its surroundings, where dust and cold gas distribute in 
a non-spherica l manner, forming a dis k-like structure of matter 
up to ~ 100/?* (IChatv & Rahouir2012h . A likely high inclination 
of the system would produce the extreme X-ray absorption and 
the eclipse-like correlation between FeKa and continuum fluxes. 

In Figure [8] we have seen that the centroid energy of FeKa 
is higher when the line is broader. When more ionized Fe goes 
along with the presence of more variety of Fe ions involved in 
the total emission, the width resulting from the blending of lines 
must depend on the centroid energy of the line, as observed. 
We have estimated the broadening produced in the lines by line 
blending by: <tb ~ E - 6.4 (keV), with E the centroid energy 
of the line in keV. We note that it is also plausible that not re¬ 
solved Fe XXV and Fe XXVI actually shift and broaden FeKa, 
producing an equivalent effect. 

More processes are also able to significantly broaden FeKa. 
We have considered Compton broadening and Doppler shifts 
as plausible candidates. Compton scattering has been proposed 
as a possible broa dening mechanism of em ission lines in neu¬ 
tron star LMXBs (iDfaz Trigo et al.l (1201 2h . for GX 13+1). For 
HMXBs, Compton broadening might also be significant, given 
the high Nh values observed (and the consequent high number 
of free electrons). However, if this process is determining the 
width of the lines, we should observe a direct correlation be¬ 
tween the absorption column and the line width. In FigureffSlwe 
can see that such a direct correlation is not p resent. Moreover, 
an inv erse correlation is plausible. Noticeably. ICackett & MIIIct 
( 1201 3h analysed three neutron star LMXBs and arrived to a simi¬ 
lar result. Therefore, Compton broadening cannot be considered 
as the main responsible for the observed width in HMXBs, al¬ 
though it is not ruled out as a modest contributor. We name ctc 
to the contribution of Compton scattering in the line broadening. 

To estimate ctc, we have used an empirica l formula accurate 
to within 30%, deriv ed from iKallmanI (Il989h and corrected in 
iBrandt & Mad d 19941) : 


o-c = 0.019 Tj-/, (1 +0.78k7’,) = 0.12 t7’/,(1 +0JSkTe) 

where Ek = 6.4 keV is the energy of FeKa, Tjh is the Thom¬ 
son optical depth and kTg the electronic temperature in keV. We 
use Trh - o-jh f ne{s) ds - ctji, Ne, where crrh is the Thomson 
cross section, tie the electron number density and the integral is 
calculated along the line of sight. Assuming solar abundances, 
a temperature kTe 1 keV and an almost completely ionized 
matter, as reasonable for galactic massive stars atmospheres, we 
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Fig. 15: Total equivalent hydrogen column (Nh) against the width of 
FeKtr. 



Transmitted 


Reflected 


Fig. 16: Simple sketch of two plausible configurations of circumstel- 
lar matter in HMXBs. In the left side. X-rays are transmitted through a 
dense medium (e.g. the strong wind of the donor), producing high 
directly correlated with the EW of FeKcr. In the right side. X-rays are 
reflected in an accretion disk producing fluorescence, and also transmit¬ 
ted through a more diffuse medium. In this case Nh is not necessarily 
correlated with the EW of EeKor. 


obtain (see e.g. equation 3.61 in lNovotn\^(ll973l) ): 


N, = 


f 


p(s)ds X(l+X) 

^:^(1+X)= Nh^ 0.7 Nh 

^ (Tc ~ 0.5 Nh{^ 0^^ cm^)eV 


Regarding Doppler shifts, they must be taken into account, 
since a velocity of more than 500 km/s (a very feasible speed, 
either in the wind either in the accretion flow), would broaden 
the lines by more than 10 eV. We name (To to any broadening 
produced by Doppler effects. 

Line blending, Compton scattering and Doppler shifts pro¬ 
duce a resultant width of: 


o-total = + o-c + *^0 

We have adopted ctd - 20,40 (eV) corresponding to velocities 
of y aj 1000,2000 km/s, which are very plausible either in the 
wind of the supergiant or either in the accretion flow. We have 
overplotted the corresponding values of crtotai in Figure^ For ev¬ 
ery observation (black squares) we have computed the expected 
value (red with ctd - 20 (eV) and green using cro - 40 (eV)). A 
vast majority of the line widths can be described in this way. 

In IGR J16318-4848, the high absorption measured (above 
2 X 10^"^ cm^) and the consequent expected Compton broaden¬ 
ing of more than 100 eV, is not congruent with the measured 
width of ~ 35 eV. This is another indication that the absorbing 
matter in this system is cold and not ionized, as already stated by 
IChatv & Rahouil(l2012h . Hence, the employed expression for de¬ 
scribing (Tc cannot apply here, since there are not sufficient free 
electrons to broaden the line by means of Compton scattering. 

In Figure |3 we show the curve of growth of FeKa. We re¬ 
quire that Nh > 2 (intrinsic absorption rather than interstellar). 
In our sample, this criterion constrains the systems in Figure 0 
to those with supergiant donors only. We observe a direct corre¬ 
lation between Nh and the EW. This correlation highlights that 
the X-ray absorption is strongly linked to the matter that emits 
FeKcr, being produced by matter in the line of sight, where the 
X-rays are absorbed, and not in other plausible regions such as 
an accretion disk, (see a sketch in Figure fTbl) . Namely, in the 


systems included in Figure|3(all with a supergiant optical star), 
FeKa is produced from the transmission of X-rays through the 
circumstellar medium, that is, either through the strong wind 
of the supergiant donor either through any structure in the line 
of sight such as ionization or accretion wakes. The hypotheti¬ 
cal reflection of X-rays in an independent medium might pro¬ 
duce an additional amount of FeKa, as observed in the BeXB 
GRO J1008-57 by iKiihnel et al.l (1201 3l) . which is not appreci¬ 
ated in the systems shown in Figure |2 

As stated before, the region where fluorescence is emitted 
must be more extended than R > R*, and consequently the wind 
of the companion star, being illuminated by the X-ray source, 
is an obvious contributor to both the absorption and the FeKa 
emission. The orbital modulation of Nh shown in Figures [12] and 
[T4|also support such an interpretation. 

Moreover, most of the observations track the numerical sim¬ 
ulations assuming a characteristic range of F values, indicat¬ 
ing that an isotropic distribution of absorbing (and FeKa emit¬ 
ting) matter is not far from reality. However, we do not ignore 
the variability and heterogeneous properties of the HMXBs en¬ 
vironment, which is very plausible to change between differ¬ 
ent sources and even between different snapshots of the same 
source. We think that this variability is reflected in the observed 
dispersion of the curve of growth and the moderate discrepancies 
regarding the simplified view of spherically distributed neutral 
matter. 

In Figure [T0| we have compared the observed values of Nh 
for SGXBs and SFXTs. We have observed that SGXBs are in 
general more absorbed sources than SFXTs. It implies that, in 
SGXBs, either the compact object orbits a denser region of the 
donor wind, or either the interaction compact object - wind mod¬ 
ifies the environment producing an enhancement of density in its 
surroundings. 

Having a look to the o rbital parameters of SFXTs (see e.g. 
Table 2 in iRomano et al.l (l2014l) ). their orbital periods lie in a 
wide range of values, from ~3 days for IGR J16479-4514, up to 
164 days for IGR J11215-5952. Some of them show high eccen¬ 
tricity. However, currently there is not a complete description 
of the orbital parameters in SFXTs. Therefore, we cannot rule 
out the possibility that, in this sample, SFXTs are less absorbed 


Article number, page 12 ofl30l 




























A. Gimenez-Garcia: Title 


than SGXBs because of the distance of the compact object to the 
donor star. In this regard, further studies of orbital parameters of 
SFXTs will be useful. 

Regarding the interaction compact object - wind, hydrody¬ 
namic simulations show that the gravitational potential of the 
compact object, and the X-ray rad iation field, can significantly 
modify the observed va lue of Nh (iManousakis & Waited 1201 it 
iManousakis et alJ2012h . In SGXBs, where the X-ray emission is 
more persistent, these effects might be stronger than in SFXTs, 
notably increasing the absorption. 

5. Conclusions 

We have performed the spectral analysis of the whole sample of 
publicly available EPN XMM-Newton observations of HMXBs 
until August, 2013, in order to describe its FeKa emission. In 
total, the study involves 46 HMXBs, 21 of them showing signif¬ 
icant FeKa emission. As expected, we dealed with a very het¬ 
erogenous set of objects and states of the sources, that must be 
properly organized. We have classified the systems in the follow¬ 
ing groups: BeXBs, SGXBs, SFXTs, y Cass analogs, HMGBs 
and peculiar sources. Furthermore, we have divided the observa¬ 
tions depending on the source behaviour in the following states: 
quiescence, flare, eclipse ingress/egress, and eclipse. With these 
criteria, we finally had a set of 108 spectra for our analysis, that 
produced the following conclusions: 

- The spectral atlas gives a qualitative description of the 
different groups of HMXBs, specially recognizable for 
SGXBs (fluorescence but not recombination Fe lines), and 
y Cass analogs (modeled by mekal models and present¬ 
ing fluorescence and recombination Fe lines). FeKtr is very 
likely an ubiquitous feature in HMXBs, but its detection 
strongly depends on the quality of the observations. SGXBs 
and SFXTs, which show the higher Nh among the HMXBs, 
tend to exhibit a more prominent Fe fluorescence. 

- The value of the centroid energy of FeKo- constrains the 
ionization state of the reprocessing material to be below 
Fe XVIII. 

- The FeKcr and continuum fluxes are well correlated, as ex¬ 
pected for the fluorescence emission of matter illuminated 
from an X-ray source. The different coefficients of correla¬ 
tion for eclipse and out-of-eclipse observations is in agree¬ 
ment with previous eclipse observations of HMXBs, in the 
sense of showing that the FeKa is produced in a region that 
ranges from the vicinity of the X-ray source to distances of 
the same order or greater than the stellar radius. 

- We confirm an inverse correlation between the X-ray lu¬ 
minosity and the EW of EeKa X-rays Baldwin effect, 
y Cass analogs do not follow this correlation. This suggests 
that the Ee Ka reprocessing scenario is fundamentally differ¬ 
ent in SGXBs and in y Cass analogs. 

- The width of EeKa is predominantly below 0.15 keV, and 
can be widely explained appealing to line blending, Compton 
broadening and moderate Doppler shifts (~1000 km/s). 

- The curve of growth in SGXBs shows a clear correlation 
of EeKa EW and Nh, indicating a strong link between the 
absorbing and the fluorescent matter. Erom numerical simu¬ 
lations, the assumption of spherically distributed absorbing 
matter is roughly correct for most of the SGXBs. 

- The Nh values observed in SGXBs are higher than in SEXTs. 
The disparity is hardly produced by chance, as shown by 
a permutation test of the sample, denoting a fundamental 
physical reason beneath. Systematic differences in the orbital 


parameters or different interaction compact object - stellar 
wind are plausible candidates to explain such a discrepancy. 
- The orbital modulation of Nh in IGR J16320-4751 and 
4U 1700-37, together with the aforementioned results, points 
to the stellar wind as the main contributor to both contin¬ 
uum absorption and EeKa emission in the case of supergiant 
donors. 

In summary, we present the most comprehensive study of 
EeKa in HMX Bs to date, complemen ting previous surveys at 
high resolution (iTorreion et alJl^lObl) . We have increased sig¬ 
nificantly the number of sources and extended the study to all 
major classes of massive binaries, leading to the already stated 
conclusions. 
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Appendix A: Tables 


Class 


Source 

Distance 

(kpc) 

lA 0535+26 

2.00 + 0.70 

2S 1845-024 

10.00 

X Persei 

0.80 + 0.14 

AXJ1820.5-1434 

8.20 + 3.50 

RXJ0146.9+6121 

1.95 

RXJ0440.9+4431 

3.30 + 0.50 

RXJ1037.5-5647 

5.00 

SAXJ2103.5+4545 

3.20 + 0.80 

Swift J045106.8-694803 

50.60 + 2.10 

V0332+53 

7.50+1.50 

4U 1538-522 

6.40+1.00 

4U 1700-37 

2.12 + 0.34 

4U 1907+09 

5.00 

EXO1722-363 

8.00 + 2.50 

GX 301-2 

3.04 

IGRJ16207-5129 

6.10 

IGR J16318-4848 

3.60 + 2.60 

IGR 116320-4751 

3.50 

IGR 116465-4507 

12.50 

SAX 11802.7-2017 

12.40 

Vela X-1 

1.90 + 0.20 

XTE10421+560 

5.00 + 4.00 


Reference 

(class) 


Reference 

(distance) 


BeXB 


^henavrin_et^ (201^ 
Bodagheeet^ (20071 
Liuet^i 

Liuet^ 

Liuet^ 

Liuet^i 

Liuet^ 

Liue^^ 
Beardmore et al. (20091 
Liu et al. (20061 



■Steele etal. (19981 
Grimme^^ (200^ 
Megie^et^ (2009) 
Kinugasa et al. (19981 
IGiarchenko,et^ 12005! 

Reig et al. (2005! 
Grimme^^ (200^ 
^a^kale^^ (2002) 
^^tlettet^ (2013) 
Neguemelaet^ (1999) 


SGXB 


IJuet^(200^ 
Liuet^(200^ 
Liuet^(200^ 
Liuet^(200^ 
Liuet^(200^ 
Nes£oli,e^^ (200^ 
Filliatre & Chatv (2004! 
Coleiro.et^ (201^ 
Wyter_&ZuritaHer^ (2007) 
Ton£i6n_e^^ (201^^ 
Liuet^(200^ 

Liu et al. (2006! 


AX1184L0-0536 

7.80 + 0.74 

IGR 100370+6122 

3.00 

IGR 111215-5952 

6.20 

IGR 116328-4726 

6.50 + 3.50 

IGR 116418-4532 

13.00 

IGRl 16479-4514 

7.50 + 2.50 

XTE 11739-302 

2.30 + 0.60 

IGR 117544-2619 

3.20+1.00 

IGR 118450-0435 

3.60 

IGR 118483-0311 

3.50 + 0.50 


Re^nolds_et^ ( 199 ^ 
Megier et al. (2009! 
Coxet^ 72005! 
Mason_e^^ (2009) 
ICager^e^^ (200^ 
Nes2oli_et^ (200^ 
Filliatre & Chatv (2004! 
Rahour_et^ (200^ 
Smith (2004! 
Tonejon_e^^ (201^^ 
^ad^ane,et^ (1985! 


SFXT 


Rmnano_et^ (20n^ 
^onzalez;Galan,et^ (201^ 
Sguera et al. (2006! 
Fiocchi_et^ (2013) 
^gueraet^ (200^ 
Roinano_et^ (2013a) 
^mith_e^^ (200^ 
Pelhzz^e^^ (200^ 
Sguera et al. (2007! 
Rahoui&Cha^ (200^ 


Nes2oli_et^ (2007) 
Neguemela^Rei^ (200£) 
Ma^etti_et^ (200^ 
Fiocchiet^ (2013) 
Drave.et^ (201^ 
Chatv et al. (2008! 
Negueruelu_et^ (200^ 
Pelhzz^e^^ (200^ 
Coe et al. (1996! 
Rahoui&Chat^ (200^ 


y-Cass-like 


yCassiopeiae 032 + 0.01 Lo^es^de^Oliveirae^^ (201^ 

HD 110432 0.39 + 0.05 Smith & Balona (2006! 

HD 119682 1.11 Rakowski et al. (2006! 

HD 157832 0.53 Lones de Oliveira & Motch (2011! 

HD 161103 1.50 + 0.50 Lones de Oliveira et al. (2006! 

HD 45314 1.50 + 0.26 Rauw etal. (2013! 

SAG 49725 2.20 + 0.60 Lo^es^de^Oliveirae^^ (201^ 

SS397 1.50 + 0.20 Lones de Oliveira (2007! 


Megieret^ (200^ 
^egieregak (2009) 
IGiMuhenko,egak (2005! 
Loges_de,01iveira^Motch (201^ 
I;Oges_de_01iveirae^^ (200^ 
Megieretak (200^ 
^oges,de_01iveiraegaL (200^ 
Riguehne_et^ (201^ 


HMGB 

LS 1+61 303 

LS 5039 

2.50 

2.50 + 0.10 

Dubus (2013! 

Dubus (2013! 

Grimm et al. (2002! 

Casares et al. (2005! 

09.5V+NS 

4U 2206+54 

2.90 + 0.20 

Blav et al. (2006! 

Riquelme et al. (2012! 

disk-fed 

Cen X-3 

10.00+1.00 

Liu et al. (2006! 

Hutchings et al. (1979! 

BH-SGXB 

Cygnus X-1 

2.10 + 0.25 

Liu et al. (2006! 

Ziolkowski (2005! 

Bl-2 

AX 11749.1-2733 

16.00 + 3.50 

Karasev et al. (2010! 

Karasev et al. (2010! 


Table A.l: Table of sources included in the sample of HMXBs. 











































































































Table A.2: Parameters of the continuum. 


Source 

Obs ID 

State 

Model 

Gain 


Po 

lA 0535+26 

0674180101 

Q 

TNi 

- 

0.87 

8.70E-01 

2S 1845-024 

0302970801 

Q 

Ts 

- 

0.73 

7.67E-01 

X Persei 

0151380101 

Q 

Ti3 

- 

1.34 

3.68E-03 

X Persei 

0151380101 

Q 

TNi 

- 

1.41 

7.44E-04 

X Persei 

0600980101 

Q 

TNi 

- 

1.28 

1.06E-02 

AX J1820.5-1434 

0511010101 

Q 

Ts 

- 

0.53 

9.92E-01 

AX J1820.5-1434 

0511010101 

Q 

Ni 

- 

0.60 

9.77E-01 

RX J0146.9+6121 

0201160101 

Q 

TNi 

- 

1.31 

8.58E-03 

RX J0440.9+4431 

0653660101 

Q 

Te 

- 

1.16 

8.55E-02 

RX J0440.9+4431 

0653660101 

Q 

Tio 

- 

1.21 

4.14E-02 

RX J0440.9+4431 

0653660101 

Q 

TNi 

- 

1.24 

2.54E-02 

RX J1037.5-5647 

0550560101 

Q 

T 9 

- 

1.02 

4.16E-01 

RX J1037.5-5647 

0550560101 

Q 

TNi 

- 

1.00 

4.94E-01 

SAXJ2103.5+4545 

0149550401 

Q 

Tio 

- 

1.26 

1.77E-02 

SAXJ2103.5+4545 

0149550401 

Q 

TNi 

- 

1.24 

2.69E-02 

Swift J045106.8-694803 

0679381401 

Q 

Ni 

- 

1.06 

3.28E-01 

V0332+53 

0506190101 

Q 

Ts 

- 

0.60 

7.58E-01 

4U 1538-522 

0152780201 

EC 

Ni 

- 

1.33 

7.12E-03 

4U 1538-522 

0152780201 

I/E 

N 2 

- 

1.44 

5.65E-04 

4U 1700-37 

0083280101 

D 

N 2 

- 

1.38 

9.69E-03 

4U 1700-37 

0083280101 

E 

Ni 

- 

0.97 

5.73E-01 

4U 1700-37 

0083280101 

Q 

Ni 

- 

1.16 

1.04E-01 

4U 1700-37 

0083280201 

E 

N 2 

- 

1.28 

1.23E-02 

4U 1700-37 

0083280201 

Q 

TN 2 

- 

0.90 

7.98E-01 

4U 1700-37 

0083280201 

Q 

N 2 

- 

1.00 

4.69E-01 

4U 1700-37 

0083280301 

Q 

N 2 

- 

1.14 

1.20E-01 

4U 1700-37 

0083280401 

EC 

Ni 

- 

1.22 

1.30E-01 

4U 1700-37 

0083280401 

I/E 

N 2 

- 

1.13 

2.01E-01 

4U 1700-37 

0600950101 

EC 

Ni 

- 

1.18 

6.97E-02 

4U 1700-37 

0600950101 

EC 

N 2 

- 

1.04 

3.58E-01 

4U 1907+09 

0555410101 

E 

TNi 

- 

1.10 

1.89E-01 

4U 1907+09 

0555410101 

Q 

TN 2 

- 

0.99 

5.16E-01 

4U 1907+09 

0555410101 

Q 

N 2 

- 

1.13 

1.30E-01 

EXO1722-363 

0206380401 

Q 

Ts 

- 

1.33 

9.84E-03 

EXO1722-363 

0206380401 

Q 

Ni 

- 

1.28 

2.30E-02 

EXO 1722-363 

0405640201 

EC 

Ts 

- 

0.71 

7.47E-01 

EXO 1722-363 

0405640201 

EC 

Ni 

- 

0.85 

5.79E-01 

EXO 1722-363 

0405640301 

Q 

Ts 

- 

0.91 

7.42E-01 

EXO 1722-363 

0405640301 

Q 

Ni 

- 

1.04 

3.67E-01 


Tl-lOicV 

{Kfi^ergjs) 


ll-03!fi39 
15.93!«4 
77 47+27.82 
72 . 49 + 27.81 


86.94^ 

12.05 


E 


J.f? 


12.37 
7.65_q j 2 
65-96!g:|J 

66 .oo ;|'3 

66.06+2|26 

■1.58*1 

11.53+0-38 
805.54:*-||^“ 
806.39+IP' 
860.26!! 
0-19!ffi 


20 04+2-24 

659 55+®''25 

JO, gg 

161.58+^^'” 


1790 09+'’42-3i 
^ '^'2''2^-545.05 
624 17+228.25 

04.^+.! /_J90 74 

1952 17+003.31 

iyjz..L /_5g9 18 

619 13+224.34 

J_igg g 2 

623.31^5 
718.06!«|« 
24.12+10-00 


313.50! 

16.64; 

16.77 


Ml 

1288.3 i+f^® 
721.87!2i 
718.59 
522.62 
528.75! 

2 77+’ 

2.82!j;| 

474 . 87 ! 3 ^ 2.48 

479.86+ 


-3>p4 

+38^.71 


.05 


•257.24 


ZNh 

{lOMmM 

r 

A 79+0.06 

^•^^-0.07 

1.40 


- 

A CO+0.08 
^•^^-0.06 

A 9 A+O.O 9 
^•■^^- 0.11 

1.63 

1.34 

7 oc+2.64 

1 78 

14 21+3-^'^ 
iH-.zi _2 03 

1.34 

0.90 

1.61 

0.78+°“ 


-Q.Q7 

1 1 7+0.05 

^•^^-0.13 

0.95 

1.04 

- 

1.04 

1.09 

A 7Q+0.06 

+'•'^-0.04 

- 

0.83!°;!J 

0.96 

0-29!°:2^ 

1.07 

0.83 

- 


kT, 

(keV) 


+0.06 

-0.08 


;8:li 

-0.14 

pO.35 

- 0.22 

fO.lO 

-0.08 


+ 0.12 

-0.13 

+0.38 

- 0.12 

+0.24 

- 0.22 

+0.13 

- 0.10 


1 94+0.05 
-0.06 
A 79 + 10.76 
^■^"^-0.92 
1 67+0.05 
^ -Q.Q7 

1 cq+0.04 

o:8:S 

^•^*^-0.23 

1 94+0.08 

0.91!; 
1.30!| 
1.38!! 
4.42+' 
1.55+' 
2.46 

2 15 

Z..1J njl 


0 40+°-°^ 
'+•^'4-0.07 


8:!5 


0 97+0-03 
-0.03 

80.48!}f[® 
11 31+^-^^ 
10.67!°'46 
15.49 
16.17!“-54 
6.02!^;® 

^•^^-a72 

23 76+0-65 

+-7-'ri_Qgg 
11.66 
97.57 
11 36+*’'^ 

1 l.JU „ jg 

15 96+0 *4 

6 39+0-52 

9 10+'“ 

17 65+'^ ® 
4 '•+’- 2 - 0.68 

24 14+*-+ 
-3.72 

13-10!g 

18.55!™ 


0.37!0'“ 
0 . 86 ! 
1.20! 
i.iO!o 
1.22 

1 1 9+0.10 
^■^^-0.03 
1 9 C+O .54 
^■'^'^-0.25 

1 03+"-“ 

4 •+'■5-0.04 

1 00+'’'“ 
7 ■+'+'-0.05 

0.51 

0.58 

0.70 

0.70 

0 75+008 
u./j 


7 cq+0.14 

^•57'-0.26 


0.45! 

1.14+: 


02 


17.78+ 


0.78!«;«^ 


-0 

^■^^-0.28 


rv 77+0.09 
' -0.08 


0.87!°;^ 

^•^^-0.03 


9 7Q+0.10 
^•/O-0.09 


4.94 


971 +0.11 
^-0.09 


kT2 kT-i 

(keV) (keV) 


10.87 


+5.49 

-3.41 


2.28 

5.82 


+0.23 

;8:i? 

-0.35 


1 c 7+0.08 
4-55-0.15 










Table A.2; Parameters of the continuum (continued). 


EXO1722-363 
EXO1722-363 
EXO 1722-363 
EXO 1722-363 
GX 301-2 
GX 301-2 
IGRJ16207-5129 
1GRJ16318-4848 
1GRJ16318-4848 
IGRJ 16318-4848 
IGRJ 16318-4848 
IGRJ 16318-4848 
IGR116318-4848 
1GRJ16318-4848 
1GRJ16318-4848 
IGRJ 16320-4751 
IGRJ 16320-4751 
IGRJ 16320-4751 
IGRJ 16320-4751 
IGRJ 16320-4751 
IGRJ 16320-4751 
IGRJ 16320-4751 
IGR 116320-4751 
IGR 116320-4751 
IGR 116320-4751 
IGR 116320-4751 
IGR 116320-4751 
IGR 116320-4751 
IGR 116320-4751 
IGR 116465-4507 
SAX 11802.7-2017 
Vela X-1 
Vela X-1 
XTEJ0421H-560 

AXJ1841.0-0536 
IGR J00370H-6122 
IGR J00370H-6122 
IGR J00370H-6122 
IGR J00370H-6122 
IGR 111215-5952 
IGR 111215-5952 


0405640401 

Q 

^5 

- 

1.09 

2.69E-01 

0405640401 

Q 

Ni 

- 

1.00 

4.77E-01 

0405640901 

Q 

Ni 

- 

1.14 

1.50E-01 

0405640801 

Q 

Ti 

- 

1.15 

1.44E-01 

0555200401 

E 

Ni 

slope 

1.29 

1.21E-02 

0555200401 

Q 

Ni 

slope 

1.82 

9.70E-09 

0402920201 

Q 

Ni 

- 

1.17 

8.16E-02 

0154750401 

Q 

Ts 

- 

1.12 

2.31E-01 

0154750401 

Q 

Ni 

- 

1.05 

3.63E-01 

0201000201 

Q 

Ts 

- 

0.89 

7.12E-01 

0201000201 

Q 

Ni 

- 

0.92 

6.36E-01 

0201000301 

Q 

Ts 

- 

0.73 

6.83E-01 

0201000301 

Q 

Ni 

- 

0.76 

6.55E-01 

0201000401 

Q 

T 5 

- 

1.12 

2.57E-01 

0201000401 

Q 

Ni 

- 

1.16 

2.12E-01 

0128531101 

Q 

Ts 

- 

0.99 

4.94E-01 

0128531101 

Q 

Tn 

- 

1.03 

4.23E-01 

0128531101 

Q 

Ni 

- 

1.03 

4.16E-01 

0201700301 

E 

Ni 

- 

0.89 

8.13E-01 

0201700301 

Q 

Ni 

- 

0.93 

6.96E-01 

0556140101 

Q 

Ni 

- 

0.99 

5.15E-01 

0556140201 

Q 

Ni 

- 

1.05 

3.47E-01 

0556140301 

Q 

Ni 

- 

1.09 

2.28E-01 

0556140401 

Q 

Ni 

- 

1.21 

6.25E-02 

0556140501 

Q 

Ni 

- 

0.84 

8.57E-01 

0556140601 

Q 

Ni 

- 

1.04 

3.73E-01 

0556140701 

Q 

Ni 

- 

1.24 

5.34E-02 

0556140801 

Q 

Ni 

- 

1.15 

1.31E-01 

0556141001 

Q 

Ni 

- 

0.93 

7.05E-01 

0164561001 

Q 

Ni 

- 

0.58 

9.23E-01 

0206380601 

Q 

Ni 

- 

1.28 

3.18E-02 

0111030101 

E 

Ni 

- 

1.35 

4.24E-03 

0111030101 

Q 

Ni 

- 

0.87 

8.65E-01 

0139760101 

Q 

Ni 

- 

0.70 

9.51E-01 


0604820301 

E 

Ni 

- 

1.39 

1.54E-03 

0501450101 

Q 

Tio 

- 

0.84 

9.23E-01 

0501450101 

Q 

Te 

- 

0.91 

7.71E-01 

0501450101 

Q 

T 9 

- 

0.97 

5.92E-01 

0501450101 

Q 

TNi 

- 

0.84 

9.24E-01 

0405181901 

E 

TNi 

- 

1.04 

3.64E-01 

0405181901 

Q 

TNi 

- 

0.98 

5.39E-01 


181.41+^^2^ 

182 

203 

287.66!215i2 
3412.87:^ 32 
1334.78!go 

94.98+2” 

II 11 + 22:83 

III c+22.66 
^^■^'^-10.31 

11 56+®*’ 

ll.JU_ig70 
11 62+2**’ 

3 61+*-2® 

n Z'^ + \6.b2 

7 56+®^ 
78.37“ 
7Q S7+™ ’’ 


918.65_ 
417.44+1^ 
474 86+**2 
920.37+|p! 
890.6321 

547.92+4*6 
443.72+*“ 


1145.80f«g’ 

:)+333.69 




418.19^ 
533.78_ 
666.20;|f27 

74.04!6-n 

624.71!|*6 

S47 4^+126.06 
-^^-112.92 

323.43+’+'® 


3.99 


-3.84 


501 72+'04.4i 
JUl./z_g4 20 

47.94+0“ 


47.70+; 

47.95+; 

48.01+' 

639.50 


76.78 


-0.97 


36.68!|28 

44.56!*;* 
30.84!*;*? 
31.45!Vn 

053 


2OO.37+0 


228.22- 


0.57 


34.91!i*o 
35.37!|| 
20.53+ 
27.23_, 
45.79+' 
32.21 


Ml 

'■M 


32.11 

38.80; 

37.64; 

33.80 

71.96 


34 90+'-67 
-1-24 

30.68!!'' 
244c 


-1 
1 +1 
- 1.01 
+ 1.58 
0.92 
1.67 
-1.26 
+3.36 


72.18+ 


2+24.40 

15 72+^®^®’ 
34.63!'“ 
39.68!0'*2 
36.88!]!* 


0 63+0-'4 

'+■'2-’-0.09 
0 97+018 
-0.07 


0.81 

1.04 


3 14 + 0.21 
-’• + ^-0.18 


A AO+0.00 
^•^■^-0.07 


50.75 


0:12 


1.07 

1.17 


1 41 +0.45 

0.45!'’ ■ 
0.65+; 
0.73! 
0.39! 
0.43! 
0.72! 
0.67! 
0.32! 
0.17! 
0.53! 
0.52! 

0 42+' 

1 - 66 ! ;J 

1-65!; 

0.57!»;| 


>.08 

>.55 


2 

^•^^-0.09 


16.27 

1.27 

- 

- 

797 1 7+19.61 
'-13.52 

- 

2 39 + 0.43 
^•o^_0.32 

- 

232.66!“4f 

206 14+I5.® 
914 99+1 l.OO 

309.57+^®*’ 
317 51* 

0.91!°;[® 

- 

- 

- 

-0.27 

- 

1.08!°;’° 

1 66 +°“ 

A 1 c+2.90 
^•^'^-0.46 

A cc+0.45 
^•'^'^-0.31 

- 

191.92+^3®' 

1 iq+ 0.45 
+ ■ + ^- 0.20 

- 

- 

24.97!*;*° 

- 

A iq+0.33 
^•^^-0.23 

- 


41.90!*-2’ 

0.87!°;“ 

0.65 

1 1 a+0.05 

Oo78:?? 
+ •'+' 2 - 0.1 3 

0 89+° ® 
-Q. 0 § 

1 02 +°’' 

+ •'+^- 0.12 


1 'yn+O.Ql 
+ ' -0.05 


1 16+°'® 
n j. 


1.04! 

0.63+ 


1 16+°°® 
+ -6^-0.08 

0 76+°°’ 
-Q.Q5 

3.08!°°® 
1.28+' 
1.87 


4 '70+0.69 
^•4'4'-Q.55 
1 00 - 1 - 0.18 
6 -°^- 0.10 


09 

; 8 :i 


1.22 


0.23 

+ 0.14 

- 0.06 







Table A.2; Parameters of the continuum (continued). 


IGRJ16328-4726 

0654190201 

Q 

Ts 

- 

1.04 

3.77E-01 

IGRJ16328-4726 

0654190201 

Q 

T9 

- 

1.03 

3.93E-01 

IGRJ 16328-4726 

0654190201 

Q 

Ni 

- 

1.07 

2.98E-01 

IGR J16418-4532 

0206380301 

Q 

Ts 

- 

0.90 

7.30E-01 

1GRJ16418-4532 

0206380301 

Q 

Ni 

- 

0.91 

7.01E-01 

IGRJ16418-4532 

0405180501 

Q 

Ts 

- 

1.11 

1.99E-01 

IGR J16418-4532 

0405180501 

Q 

Ni 

- 

1.10 

2.26E-01 

IGRJ16479-4514 

0512180101 

EC 

Ni 

- 

0.95 

5.92E-01 

IGRJ16479-4514 

0512180101 

I/E 

N 3 

- 

0.74 

9.67E-01 

XTEJ1739-302 

0554720101 

E 

Ts 

- 

1.25 

5.66E-02 

XTEJ1739-302 

0554720101 

E 

Tu 

- 

1.08 

2.88E-01 

XTE J1739-302 

0554720101 

E 

Ni 

- 

1.20 

9.43E-02 

XTEJ1739-302 

0554720101 

Q 

Ts 

- 

1.20 

2.51E-01 

XTE J1739-302 

0554720101 

Q 

Tn 

- 

1.39 

1.26E-01 

XTE J1739-302 

0554720101 

Q 

Ni 

- 

1.52 

7.23E-02 

XTE J1739-302 

0561580101 

E 

Ts 

- 

0.96 

6.04E-01 

XTEJ 1739-302 

0561580101 

E 

Tu 

- 

0.84 

8.79E-01 

XTEJ1739-302 

0561580101 

E 

Ni 

- 

0.88 

8.05E-01 

XTEJ 1739-302 

0561580101 

Q 

Ts 

- 

1.24 

6.93E-02 

XTEJ 1739-302 

0561580101 

Q 

Tu 

- 

1.02 

4.33E-01 

XTEJ 1739-302 

0561580101 

Q 

Ni 

- 

1.09 

2.77E-01 

IGRJ 17544-2619 

0148090501 

E 

Ts 

- 

0.90 

7.16E-01 

IGRJ 17544-2619 

0148090501 

E 

Tu 

- 

0.81 

8.91E-01 

IGR J17544-2619 

0148090501 

E 

Ni 

- 

0.87 

7.83E-01 

IGRJ 18450-0435 

0306170401 

E 

TNi 

- 

1.04 

3.51E-01 

IGR J18450-0435 

0306170401 

Q 

TNi 

- 

1.05 

3.39E-01 

IGRJ 18483-0311 

0406140201 

Q 

Ts 

- 

1.17 

1.98E-01 

IGRJ 18483-0311 

0406140201 

Q 

T 9 

- 

1.28 

9.12E-02 

IGR J18483-0311 

0406140201 

Q 

Ni 

- 

1.39 

4.16E-02 


99 A 7 + 45.82 
jz..u/_25.42 

16.67![;| 

- 

2.14^ 

32.35!4“o 

22 . 21 !';®« 

25.13!“® 

- 

4.48^ 

99 94+4643 

1 41 +0.18 
^■^^-0.13 

- 

122.35!«« 

10.73![;« 

- 

2 . 1 E 

124.60!“3 

78.92“-| 

82.17““ 

17.64!f| 

3.92+°® 

1 99 + 0.22 
^■■^'^-0.15 

IJT 

0 oc+S:?^ 

tS >49 f\ c-j 

1.49!°® 

_ 

5 

'^•■^^- 3.10 


1 75+^-^^ 
^■^^-0.26 

- 


0.82 

0.12 

M 


978 ■ 1 ' 7 + 225.37 

1 71^”'' 

0.16^ 

A 1 q+0.13 
U.l»_g.o9 

5.68!|® 

5.96!^! 

6.00!^;*^ 

1 14+0.74 

16:583f 

131.25!|'4' 

1.25+'!™ 


1.26!; 

1.29+: 


69.51 
1 51+°-^* 

38 

■^•^^-0.32 
A 90+0.46 
^•^^-0.33 

;157?:§A 

5.08+!^® 
26 


8 :i 

0.42 

+0.50 

-0.41 

+0.29 

0.25 

+0.38 


2.24+' 
5.12 
5.31 
2.28 

4 94 

5 52^8:?? 
''•'’^-0.42 

^•^'+-0.28 

4 58 +°+^° 

^•'’<’-0.43 

5 17+0-75 

-0.48 

2.18!“! 
2 47 + 0.58 
^•^'-0.27 

4 64+' " 

^•"^- 0.75 
7 47+1-28 

11 40 + 0-53 

i I.HU 0 33 


1.12 


1 SI+0-13 
^■°7_o.lO 


9 /lQ+0.46 
^■^^-0.30 


1 49+0.11 
^■^*^-0.09 


1 00+015 
^■^^-0.10 


1 Q9+0.21 
^■^^-Q.14 

I 

19 


0.84;^; 

0.79^; 


9 44+0.15 
^■^^-0.14 


+ 0.11 

0.10 

+0.45 

0.59 

+0.15 

0.13 

+0.07 

-0.07 


1 97+0.05 
^ ^ -0^05 

10 

g4 


1 

^•^^- 0.12 
9 Q9+2.55 
*^•^^-1.37 

1 57+0.06 

.4274^8 

-8.96 


1 49+0.06 
^•^^-0.05 
10.85“;! 


1 99+0.07 
8.60!?;| 


1.58!“® 

4+0!23 


L74_^ 

L36!j 

2.07+' 


10 

25 

0.22 


y Cassiopeiae 

0201220101 

Q 

T 3 

slope 

1.39 

1.42E-03 

y Cassiopeiae 

0651670201 

Q 

T 4 

- 

1.22 

3.62E-02 

y Cassiopeiae 

0651670201 

Q 

T 3 

- 

0.89 

8.39E-01 

y Cassiopeiae 

0651670301 

Q 

T 3 

- 

1.43 

4.40E-04 

y Cassiopeiae 

0651670401 

Q 

T 3 

offset 

1.15 

1.09E-01 

y Cassiopeiae 

0651670501 

Q 

T 3 

offset 

1.18 

6.73E-02 

HD 110432 

0504730101 

Q 

Ni 

- 

1.24 

2.81E-02 

HD 110432 

0504730101 

Q 

TNs 

- 

1.03 

3.78E-01 

HD 119682 

0551000201 

Q 

Ti 

- 

1.11 

2.06E-01 

HD 119682 

0551000201 

Q 

TN 3 

- 

1.18 

1.06E-01 

HD 157832 

0551020101 

Q 

Ti 

- 

1.07 

3.12E-01 

HD 161103 

0201200101 

Q 

Ti 

- 

1.20 

9.76E-02 


0.24!°;°® 

0.56 

_ 

0.96 

5.27 

33.04 

A 97+0.06 

02 + 8 : 8 ? 

+-Q5 

0.24+°°® 

0 23+°-°' 

0 42+0-05 

'+•^^-0.04 

0.56 

“ 

28.62!“4 

0.74 

0.64 

3.42 

3.95 

31.13 

40.83 

0.33+11 

0.56 


0.66 

3.58 

38.62 

0.23^1 

0.56 

_ 

0.66 

3.40 

38.38 

0 21 + 8 - 8 ° 
02 + 8:88 

01 + 8:88 

-Q-Qi 

0-16!£ 

0 27+0.02 
-0.02 

0 11+0-01 
i^-0.02 

1.81 

1.81 

1 99+0.02 
^■'^^-0.02 

0 94+0.18 
-0.15 

1 63+°'® 

8 69+'-'® 

A 99+0.00 

^•^^-Q.Ol 

14 53+° °'' 

- 

0-08!:°8 

0.34 

0.94 

- 

6 50+° °° 

U.JU 0 73 

1.03 

7.05 

; 





Table A.2; Parameters of the continuum (continued). 


HD 45314 

SAG 49725 

0670080301 

0201200201 

Q 

Q 

T 2 

Ti 

- 

0.96 

1.31 

5.63E-01 

7.39E-02 

0 30+° ” 

0 52+° ” 

0 55+° ” 

U.JJ_g29 

0 17+007 

-QQ5 

0 17+0.06 
-QQ5 

0 20+° “ 

-QQ6 

0 20+° “ 

0 13+006 

0 14+0.06 

-QQ6 

0 19+0“ 

^-0.07 

0.73 

0.77 

SAG 49725 

0201200201 

Q 

TNs 

- 

1.19 

1.67E-01 

0.77 

SS397 

0122700101 

Q 

Ti 

- 

0.84 

7.90E-01 

1.37 

SS397 

SS397 

0122700101 

0122700201 

Q 

Q 

Ni 

Ti 

- 

1.00 

0.89 

4.67E-01 

6.88E-01 

C9 

d d 

SS397 

0122700201 

Q 

Ni 

- 

0.93 

6.04E-01 

1.37 

SS397 

0122700301 

Q 

Ti 

- 

1.03 

4.21E-01 

1.37 

SS397 

0122700301 

Q 

Ni 

- 

1.08 

3.35E-01 

1.37 

SS397 

0122700501 

Q 

Ti 

- 

1.01 

4.49E-01 

1.37 

SS397 

0122700501 

Q 

Ni 

- 

0.95 

5.45E-01 

1.37 

LS 1+61 303 

0207260101 

Q 

TNi 

- 

1.09 

2.25E-01 

6 27+° °^ 

5 89+0.13 

6 31+° ® 

6 00+° ” 

6 39+016 

^ -0.11 

11 94+0.24 

Q 07+0.18 

-Q-15 

^•^"^-0.12 

5 

''^•^^-0.52 

5 10+0.56 

5.07!;« 

5.08:°g 

5 15+°” 

'^•^'^-0.58 

6 

6 

'^•^'^-0.59 

'^•'^^-0.38 

■^•^^-0.38 

0.90 

LS 1+61 303 

0505980801 

Q 

TNi 

- 

1.03 

4.04E-01 

0.90 

LS 1+61 303 

0505980901 

Q 

TNi 

- 

1.19 

8.86E-02 

0.90 

LS 1+61 303 

0505981001 

Q 

TNi 

- 

1.09 

2.44E-01 

0.90 

LS 1+61 303 

LS 1+61 303 

0505981101 

0505981201 

Q 

Q 

TNi 

TNi 

- 

0.92 

0.99 

7.24E-01 

5.09E-01 

0.90 

0.90 

LS 1+61 303 

0505981301 

Q 

TNi 

- 

1.06 

3.03E-01 

0.90 

LS 1+61 303 

0505981401 

Q 

TNi 

- 

1.16 

L03E-01 

0.90 

LS 5039 

0151160201 

Q 

Tn 

- 

1.18 

9.13E-02 

1.03 

LS 5039 

0151160201 

Q 

Ni 

- 

0.97 

5.76E-01 

1.03 

LS 5039 

0151160301 

Q 

Tn 

- 

1.21 

6.01E-02 

1.03 

LS 5039 

0151160301 

Q 

Tio 

- 

1.20 

6.97E-02 

1.03 

LS 5039 

0151160301 

Q 

Ni 

- 

1.11 

1.99E-01 

1.03 

LS 5039 

0202950201 

Q 

Tn 

- 

1.15 

L14E-01 

1.03 

LS 5039 

0202950201 

Q 

Tio 

- 

1.10 

2.09E-01 

1.03 

LS 5039 

0202950201 

Q 

Ni 

- 

0.93 

6.99E-01 

1.03 

LS 5039 

0202950301 

Q 

Tn 

- 

1.19 

8.41E-02 

1.03 

LS 5039 

0202950301 

Q 

Ni 

- 

0.92 

7.12E-01 

1.03 

4U 2206+54 

0650640101 

Q 

TNi 

- 

1.19 

5.04E-02 

294 46 

0 03 

Cen X-3 

Cen X-3 

0111010101 

0111010101 

I/E 

I/E 

TNa 

N 2 

- 

1.28 

1.27 

L24E-02 

L39E-02 

2242 07+477.52 
^^^^■^'-430.86 
2240 

10191.40!™ 

8572.45!g°;| 
5061 89+'®'^“ 
3453.30+^++ 
5009.50!”“® 

'^•^^-0.05 

0 qc+0.74 

0 77+0-02 
-0.01 

Cygnus X-1 
Cygnus X-1 

0202401101 

0202401201 

Q 

Q 

TNa 

TN 3 

- 

1.12 

1.13 

L54E-01 

L39E-01 

Cygnus X-1 

0202760201 

Q 

Ni 

- 

1.49 

2.94E-03 

0.86 

Cygnus X-1 

0202760301 

Q 

Ni 

- 

0.96 

5.82E-01 

0.86 

Cygnus X-1 

0202760401 

Q 

Ni 

- 

1.72 

6.18E-05 

0.86 

Cygnus X-1 

0202760501 

Q 

Ni 

- 

1.17 

L38E-01 

4842 41+123^95 

0.86 

Cygnus X-1 
Cygnus X-1 

0500880201 

0500880201 

D 

D 

TN 4 

N 2 

; 

1.19 

1.19 

6.15E-02 

5.54E-02 

1899.40”®” 

1899 66+®° “ 

A C71+0.26 
^•'^■^-0.34 

■^'^^-0.45 

Cygnus X-1 

0500880201 

Q 

Ni 

- 

1.46 

3.73E-04 

2679 69+“° °° 

^"'^■'^^-600.04 

0.71 


0.82+0;40 

1 70+0.24 
^■^■^-0.16 

1 . 58!«:[2 


1.42-*-o 


1.98!° 

2.34;^; 

2.47+; 

2.28+ 

2.08 

1.89+ 

2.07 

2.06_o 


05 




8:4^ 

11 

11 

12 

11 

17 


1 60+'’°'^ 


1.48 


+0.04 

-0.04 


1 ^4+0.03 
^ ■-'+^-0.03 


1 59 + 0.04 
^ -"^-0.04 


1 01+004 
^■'4r_o.04 

0.68 

0.71 

2 57+0.06 

2 51+0.05 

-8:81 
).01 
).01 

!:84 
).01 
).02 

!: 8 ? 


2.29+ 

2.07+ 

2 . 10 + 

2.18+ 

1.43+ 

1.42 


0.03 

1.55 


1.03 
0.82 
1 07 + 0.92 

n 'mi 

^•^■^-1.14 


13.29 

18.86l+«4 


10-03!’:“ 

6-92!l:i 

28.oi!5;-8o 


2 45+1.65 
^ ^-0.26 
1.78+°” 

1.85 

I. 79 
2.17_, 
2.58+' 
2.07 
2-12_„ 

II. 87+ 


i8:i3 

78:41 

?:?t 

8:11 

?:il 

"’?6 
1.18 


18 . 50 !: 4-55 

2QJ+0.42 0 2^+0.16 

^•^^- 0.16 ^'^^- 0.07 


14.66”“ 

7 1 A+0.26 A 07+O.O8 

^*^^- 0.16 ^•^^- 0.06 

Q 24+1.38 
^•^^-0.93 


1 57+0.04 
-0.03 

0.94 


0.46 

0.43 


+0.01 

-Q.Ol 

+ 0.00 


- 0.01 


0 39-^0-02 
^•*^^-0.02 

4.63 







Table A.2; Parameters of the continuum (continued). 


Cygnus X-1 

0610000401 

Q 

TN4 

- 

1.13 

1.43E-01 

3840.30+^^’'45 
290 03"*^ 


1 

^■^^- 0.01 

0 57 + 0-02 
- 0.01 

- 

- 

AX J1749.1-2733 

0510010401 

Q 

T5 

- 

0.92 

6.80E-01 

22.81+2™ 


^•^^-0.16 

- 

- 

AX J1749.1-2733 

0510010401 

Q 

Ni 

- 

0.95 

6.06E-01 

OQ'i 9Q+156.74 

z.yj.z.y^2i29 

31.98!2| 

1 99 + 0.18 
^■■^^-0.16 

- 

- 

- 


Table A.2: Parameters of the continuum, including observation ID, state of the source, model (see Table[T]i, artificial gain (offset or slope), reducednull hypothesis probability 
(Pq), luminosity between 1-10 keV (Li^iokev), total Nh (adding every absorption component), photon index (T), and temperature of thermal components (kT,): bbody, diskbb, 
bremss, mekal or cemekl (depending on the model). Parameters frozen or unbounded are included without an error estimation. Therefore they are not used in the plots and the 
subsequent discussion. The possible states are; quiescence(Q), flare(F), eclipse ingress/egress(I/E), eclipse(E) and dip(D). 



Table A.3; FeKo- parameters. 


Narrow FeKa 


Broad Fe feature 


Source 

Obs ID 

State 

Po 

Energy 

Width 

EpeKa 

EW 

Signif. 

Energy 

Width 

EpeKa 

EW 

Signif. 





(keV) 

(eV) 

(lO^'y/s) 

(keV) 

(O') 

(keV) 

(eV) 

(lO^'y/s) 

(keV) 

(o-) 

lA 05354-26 

0674180101 

Q 

8.70E-01 

- 

- 

< 6.92E - 02 

< 0.06 

- 

- 

- 

- 

- 

- 

2S 1845-024 

0302970801 

Q 

7.67E-01 

- 

- 

< 1.22E-H00 

< 0.60 

- 

- 

- 

- 

- 

- 

X Persei 

0151380101 

Q 

3.68E-03 

- 

- 

< 2.11E-01 

< 0.03 

- 

- 

- 

- 

- 

- 

X Persei 

0600980101 

Q 

1.06E-02 

- 

- 

< 1.76E-01 

< 0.02 

- 

- 

- 

- 

- 

- 

AXJ1820.5-1434 

0511010101 

Q 

9.92E-01 

- 

- 

< 4.03E - 01 

< 0.21 

- 

- 

- 

- 

- 

- 

RX J0146.9H-6121 

0201160101 

Q 

8.58E-03 

- 

- 

< 4.25E - 02 

< 0.06 

- 

- 

- 

- 

- 

- 

RX J0440.9H-4431 

0653660101 

Q 

8.55E-02 

- 

- 

<4.73£-01 

< 0.06 

- 

- 

- 

- 

- 

- 

RX J1037.5-5647 

0550560101 

Q 

4.16E-01 

- 

- 

< 1.72E-01 

< 0.13 

- 

- 

- 

- 

- 

- 

SAX J2103.54-4545 

0149550401 

Q 

1.77E-02 

6-42:8:8^ 

54.00 

3.87£ + 00:+“ 

A Q44-0.01 

3.27 

- 

- 

- 

- 

- 

Swift J045106.8-694803 

0679381401 

Q 

3.28E-01 

- 

- 

< 1.80E-H01 

< 0.20 

- 

- 

- 

- 

- 

- 

V03324-53 

0506190101 

Q 

7.58E-01 

- 

- 

< 1.14E-01 

< 2223.78 

- 

- 

- 

- 

- 

- 

4U 1538-522 

0152780201 

EC 

7.12E-03 

6.38:®:82 

6.40:®;« 

0.00 

9.58£-0i:®“[ 
3.42£ + 00:+^®f:®® 

A 1 0 4-0.04 
^■^^-0.04 

^■^■^-0.01 

8.70 

6 54+0.02 
^■'^^-0.02 

271 

z/ i.OJ_jg 92 

3.15£ + 00:l;«ft®8 

A 09+0.14 

19.95 

4U 1538-522 

0152780201 

l/E 

5.65E-04 

24.24 

6.46 

- 

- 

- 

- 

- 

4U 1700-37 

0083280101 

D 

9.69E-03 

- 

- 

< 6.03E -H 00 

< 0.27 

- 

- 

- 

- 

- 

- 

4U 1700-37 

0083280101 

E 

5.73E-01 

6.48:®;“ 

0.02 

6.86£ + 00:«^f:®® 

^■^■^-0.02 

3.56 

- 

- 

- 

- 

- 

4U 1700-37 

0083280101 

Q 

1.04E-01 

- 

- 

< 1.23E-H01 

< 0.13 

- 

- 

- 

- 

- 

- 


4U 1700-37 
4U 1700-37 
4U 1700-37 
4U 1700-37 
4U 1700-37 
4U 1700-37 
4U 1907H-09 
4U 1907H-09 
EXO1722-363 
EXO1722-363 
EXO 1722-363 
EXO 1722-363 
EXO 1722-363 
EX01722-363 
GX 301-2 
GX 301-2 
IGRJ16207-5129 
IGRJ16318-4848 
IGRJ 16318-4848 
IGR J16318-4848 
IGR J16318-4848 
IGRJ 16320-4751 
IGRJ 16320-4751 


0083280201 

0083280201 

0083280301 

0083280401 

0083280401 

0600950101 

0555410101 

0555410101 

0206380401 

0405640201 

0405640301 

0405640401 

0405640901 

0405640801 

0555200401 

0555200401 

0402920201 

0154750401 

0201000201 

0201000301 

0201000401 

0128531101 

0201700301 


E 

Q 

Q 


1.23E-02 

7.98E-01 

1.20E-01 


EC 1.30E-01 


1/E 

EC 

E 

Q 

Q 

EC 

Q 

Q 

Q 

Q 

E 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

E 


OlE-01 

97E-02 

89E-01 

16E-01 

84E-03 

47E-01 

42E-01 

69E-01 

50E-01 

44E-01 

21E-02 

70E-09 

16E-02 

31E-01 

12E-01 

83E-01 

57E-01 

94E-01 

13E-01 


6 45+001 


6.44+; 
6.45+ 
6.49+ 
6.48+ 
6.40+ 
6.43+ 
6.43+ 
6.40+ 
6.41 + 
6.42+ 
6.41 + 
6.42+ 
6.43+ 
6.37+ 
6.38+ 
6.42+ 
6.40+ 
6.41 + 
6.42+ 
6.41 + 


1.00 

1.03 

I.Q2 

1.02 

!: 8 ? 

>.01 

>.05 


>.01 

>.00 

>.00 

>.01 

m 


1 
1 

0.01 


6 41 + 0.02 
^■^^-0.01 


90 . 60!‘'-35 

125.45!tp 

112 97+®-82 
-25.76 

0.00 

38.57!«3 

34.05 
56 22+^^°® 

ju.z,z,_4g4i 

46.25 

124.04!||® 
54.05 ' 

45 10+20.89 

17.75 

50.70!23.3v 

25.43!^f 

0.14 

0.00 

35.48:®;™ 

26.59 

32.52 

27.11 


31.55 


2ME- 


.13E+01 
-1.31£+Q1 

9.07£ + 00:JpOj> 
8.08£ + 00:|™ 
3.547? -+ nn+i.sofi+oo 


0 11+000 
3 1-0.03 


H-OO:!^ 

' H- oo:[ 
7.28E H- 00:+ 
2.47 E + 00:f 
6.04E H- OO:}' 
5.29 £h-00:’' 
1.80£h-01:|' 
1.02£h-00:^' 
1.36£h-01:'' 
1.46£h-01:|' 
2.99E + OO:: 
1.57E 
2.19E 


■mm 


: + 0i:'- 
'+02:^- 
1.41 £h- 02:4- 
4.19£-0i:f 
2.55E 


1.05E 
3.96E 
4.56E 
< 3 
1.92E 


H- 02:« 
H- 02:^ 
+02:| 
+ 01-4' 
!.78£' -H 
+ 01!! 


.45£+Q 

13E+0 

£4-0 

;£-0 

77£-01 

10£4-00 

59£+00 

56£-01 

56£-Ql 

67£4-01 

09£+01 

23£+00 

50£-Ql 

44£4-01 

nm 

lifts? 

mm 

mm 

mm 

.39£4-Q2 

.07£+02 

.88£+Q2 

.51£4-02 

.31£+01 

00 

.84£+01 

.17£4-01 


0.11!; 

0.07!; 

1.25!; 

0.13!i 

0.92!; 

0.04!; 

0.05!; 

0.15-"! 


)2 

).03 

).02 

).02 

>.27 


8:8^ 
0.04 
O.Ol 
0.01 
O.Ol 
0.01 
1).02 
0.03 

2 49 + 0.71 


0.14+ 
0.23+ 
0.05+ 
0.17+ 
0.46+ 
0.44+ 
0.03+ 
2.47+ 
1.12 ' 


:;8? 


1.04 

1.66 


1.03 
1.00 
1.02 
1.03 
I.Q2 
1.02 

1.21 
1.13 

9 7c+?;l^ 
■ =8:8 
3-3-3-0.16 
< 0.28 

0 13+®°' 
3'-3-3-0.02 


19.07 

7.25 

8.15 

12.50 

6.82 

>25.5 

8.55 

11.60 

13.19 
7.92 
8.99 
14.26 
4.02 

14.20 
>25.5 
>25.5 

2.71 

>25.5 

>25.5 

10.97 

22.68 

5.38 


6 54+0 04 
O’■-’+^-0.03 


410.9i*fsm 


I.OIE 


aa+4.84£-01 

^^-3.69£-01 


0.19!g;g2 21.33 











Table A.3: FeKo' parameters (continued). 


IGRJ16320-4751 

0201700301 

Q 

6.96E-01 

IGRJ16320-4751 

0556140101 

Q 

5.15E-01 

IGRJ 16320-4751 

0556140201 

Q 

3.47E-01 

IGRJ 16320-4751 

0556140301 

Q 

2.28E-01 

IGRJ 16320-4751 

0556140401 

Q 

6.25E-02 

IGRJ 16320-4751 

0556140501 

Q 

8.57E-01 

IGRJ 16320-4751 

0556140601 

Q 

3.73E-01 

IGRJ 16320-4751 

0556140701 

Q 

5.34E-02 

IGRJ 16320-4751 

0556140801 

Q 

1.31E-01 

IGRJ 16320-4751 

0556141001 

Q 

7.05E-01 

IGRJ16465-4507 

0164561001 

Q 

9.23E-01 

SAX JI802.7-2017 

0206380601 

Q 

3.18E-02 

Vela X-1 

0111030101 

F 

4.24E-03 

Vela X-1 

0111030101 

Q 

8.65E-01 

XTE J04214-560 

0139760101 

Q 

9.51E-01 


6.42+0 

6.44+ 

6.42+ 

6.42+ 

6.43+ 

6.41 + 

6.42+ 

6.42+ 

6.42+ 

6-42!o 


6.45 

6.45 

6.43 


+ 0.02 



-0.03 


0.00 

0.00 

26.65 

0.00 

68.31!^o.44 

0.00 

34.48!‘J« 

45 . 43!|55 

52 95+20.fe 
21.72 


55.91 

0.00 


7.80£ 

1.83£ 

2.92E 

2.23E 

1.64£ 

1.13£ 

3.12E 

2.83E 

1.67E 

1.12E 

< 5 

< 1 
1.21E 
4.00E 
2.18E 


, nr)+7.17£+00 
'^'^-4.60£+QQ 

, Ai +1.62£+0I 

^^-1.39£4-Q1 

. Ai +1.60£+01 

+ '^^-7.30£+00 
, ni+3.25£+01 
+ '4i_2 29£+Qi 
, AI +2.46£+01 
+ '^^-1.77£+Q1 

, AI +1.61£+01 

+ ^ -l.Ql£+01 
, AI +1.82£+0l 
+ '^^-1.05£+01 

.21E + 00 
.36£ + 01 

, aa+3.78£+00 
+ '^'^-2.56£+QQ 
, aa+1.50£+00 

_ AI +1.61£+00 

'^^-2.69£-01 


0 . 11 _, 
0 . 22 + 
0.19+ 
0.15+ 
0.17+ 
0.15+ 
0 . 20 + 
0.41 + 
0.19+ 
0.15+ 


+ 0.01 
81 
i? 


1.03 

1.02 

l:i 

1 : 8 ? 

1.04 

1.04 

1 : 8 ? 

1 : 8 ? 

0.02 


< 0.50 

< 0.14 


0.07 

0.06 

0.15 


+ 0.02 



-0.04 


9.58 

- 

- 

- 

- 

- 

15.57 

- 

- 

- 

- 

- 

6.60 

- 

- 

- 

- 

- 

5.52 

- 

- 

- 

- 

- 

17.17 

- 

- 

- 

- 

- 

7.71 

- 

- 

- 

- 

- 

15.70 

- 

- 

- 

- 

- 

10.46 

- 

- 

- 

- 

- 

19.00 

- 

- 

- 

- 

- 

8.66 

— 

— 

— 

— 

- 

9.75 

- 

- 

- 

- 

- 

19.86 

2.65 

6.62+°;;« 

6.56:g 

1220.66:ig™ 

402.38!£V^ 

4.29£ + 0l!j;pf:«; 

7.08E-01!™! 

0 75+°'^ 
^+-0-09 
0 43+^ ® 

>25.5 

3.51 


AXJ1841.0-0536 

0604820301 

F 

1.54E-03 


0.00 

1.67£ + 00![;[«f:0“ 

0 

2.78 

- 

- 

- 

- 

- 

IGR1003704-6122 

0501450101 

Q 

9.23E-01 

- 

- 

< 2.91E - 01 

< 0.06 

- 

- 

- 

- 

- 

- 

IGRJl 1215-5952 

0405181901 

F 

3.64E-01 


62 01+^^'^° 
Ui.ui 45 59 

4.16£ + 00![;™ 

0 

^■^'^-0.02 

3.49 

- 

- 

- 

- 

- 

IGR J11215-5952 

0405181901 

Q 

5.39E-01 

- 

- 

< 1.16F-H00 

< 0.13 

- 

- 

- 

- 

- 

- 

IGRJ 16328-4726 

0654190201 

Q 

3.77E-01 

- 

- 

< 4.94F - 01 

< 0.09 

- 

- 

- 

- 

- 

- 

IGRJ16418-4532 

0206380301 

Q 

7.30E-01 

- 

- 

< 1.34F-H00 

< 0.07 

- 

- 

- 

- 

- 

- 

IGR J16418-4532 

0405180501 

Q 

1.99E-01 

- 

- 

< 6.39F-01 

< 0.05 

- 

- 

- 

- 

- 

- 

IGRJ16479-4514 

0512180101 

EC 

5.92E-01 

^■^^-0.00 

0.00 

^^_3.34£_oi 

0 

^■^'^-0.14 

6.36 

- 

- 

- 

- 

- 

IGRJ16479-4514 

0512180101 

I/E 

9.67E-01 

6.40 

0.00 

1.72F + 00 

0.03 

1.48 

- 

- 

- 

- 

- 

XTEJ1739-302 

0554720101 

F 

5.66E-02 

- 

- 

< 3.81F-02 

< 0.21 

- 

- 

- 

- 

- 

- 

XTEJ1739-302 

0554720101 

Q 

2.51E-01 

- 

- 

< 1.83F-02 

< 2.27 

- 

- 

- 

- 

- 

- 

XTEJ1739-302 

0561580101 

F 

6.04E-01 

- 

- 

< 9.58F - 02 

< 0.12 

- 

- 

- 

- 

- 

- 

XTEJ1739-302 

0561580101 

Q 

6.93E-02 

- 

- 

< 2.35E - 02 

< 0.17 

- 

- 

- 

- 

- 

- 

IGR J17544-2619 

0148090501 

F 

7.16E-01 

- 

- 

< 2.39F-01 

< 0.12 

- 

- 

- 

- 

- 

- 

IGRJ 18450-0435 

0306170401 

F 

3.51E-01 

- 

- 

< 2.24F -H 00 

< 0.13 

- 

- 

- 

- 

- 

- 

IGRJ 18450-0435 

0306170401 

Q 

3.39E-01 

- 

- 

< 4.27F - 01 

< 0.14 

- 

- 

- 

- 

- 

- 

IGRJ18483-0311 

0406140201 

Q 

1.98E-01 

- 

- 

< 4.23F - 02 

< 0.27 

- 

- 

- 

- 

- 

- 


y Cassiopeiae 

0201220101 

Q 

1.42E-03 

6 43+0.01 

6 39+0.03 

5 44+0.25 

” -Q.Q3 

6 41+0.02 

-Q.Q2 

6 39+0.02 

6 44+0.03 
^■^^-0.04 

0.00 

y Cassiopeiae 

0651670201 

Q 

3.62E-02 

23.94 

y Cassiopeiae 

0651670301 

Q 

4.40E-04 

6.39 

y Cassiopeiae 

0651670401 

Q 

1.09E-01 

1.50 

y Cassiopeiae 

0651670501 

Q 

6.73E-02 

0.00 

HD 110432 

0504730101 

Q 

2.81E-02 

0.14 

HD 119682 

0551000201 

Q 

2.06E-01 


- 

HD 157832 

0551020101 

Q 

3.12E-01 

6.40 

0.00 

HD 161103 

0201200101 

Q 

9.76E-02 

6.40 

0.00 

HD 45314 

0670080301 

Q 

5.63E-01 

6 39+^ ’^^ 

0.00 


6.86£ - 04 
1.08£-03 
1.00£-03 
1.53£-03 
8.55£ - 04 
1.21£-03 


+2.62£-04 

-1.95£-Q4 

+5.03£-04 

imm 

mux 

mm 

mm 

-5.97£-04 


< 1.53E - 03 


1 '^9 Z7 _ 

l.DA£L ^■^_9.18£_Q4 

4.68£-03!”^f:| 
5.45£ - 03 


0.03 

0.04 

0.04 

0.05 

0.04 

0.05 


+ 0.00 

- 0.00 

+ 0.01 

; 8:?1 

- 0.01 

+ 0.01 

; 8 : 8 ? 

- 0.01 

+ 0.01 

- 0.02 


< 0.63 


A jA+o.oe 

0 09+° °^ 
' 4 '' 42 '_ 0.07 

0.17 


13.80 

7.82 

7.41 

10.26 

7.91 

5.36 

2.47 

2.18 

1.98 












Table A.3: FeKo' parameters (continued). 


SAG 49725 

0201200201 

Q 

7.39E-02 

- 

- 

< 3.26E - 02 

< 0.71 

- 

- 

- 



- 

- 


- 

SS397 

0122700101 

Q 

7.90E-01 

- 

- 

< TOOL-02 

< 0.42 

- 

- 

- 



- 

- 


- 

SS397 

0122700201 

Q 

6.88E-01 

- 

- 

< TOOL-02 

< 0.43 

- 

- 

- 



- 

- 


- 

SS397 

0122700301 

Q 

4.21E-01 

- 

- 

< l.lOL-02 

< 0.57 

- 

- 

- 



- 

- 


- 

SS397 

0122700501 

Q 

4.49E-01 

- 

- 

< 2.69E - 02 

< 1.31 

- 

- 

- 



- 

- 


- 

LS 1+61 303 

0207260101 

Q 

2.25E-01 

- 

- 

< 2.16L-02 

< 0.04 

- 

- 

- 



- 

- 


- 

LS 1+61 303 

0505980801 

Q 

4.04E-01 

- 

- 

< 5.80L - 02 

< 0.12 

- 

- 

- 



- 

- 


- 

LS 1+61 303 

0505980901 

Q 

8.86E-02 

- 

- 

< 8.59L - 02 

< 0.14 

- 

- 

- 



- 

- 


- 

LS 1+61 303 

0505981001 

Q 

2.44E-01 

- 

- 

< 9.20L - 02 

< 0.17 

- 

- 

- 



- 

- 


- 

LS 1+61 303 

0505981101 

Q 

7.24E-01 

- 

- 

< 4.12E - 02 

< 0.08 

- 

- 

- 



- 

- 


- 

LS 1+61 303 

0505981201 

Q 

5.09E-01 

- 

- 

< 1.16L-01 

< 0.10 

- 

- 

- 



- 

- 


- 

LS 1+61 303 

0505981301 

Q 

3.03E-01 

- 

- 

< 5.54L - 02 

< 0.06 

- 

- 

- 



- 

- 


- 

LS 1+61 303 

0505981401 

Q 

1.03E-01 

- 

- 

< 6.98L - 02 

< 0.12 

- 

- 

- 



- 

- 


- 

LS 5039 

0151160201 

Q 

9.13E-02 

- 

- 

< 4.72E - 02 

< 0.10 

- 

- 

- 



- 

- 


- 

LS 5039 

0151160301 

Q 

6.01E-02 

- 

- 

< 1.06L-01 

< 0.22 

- 

- 

- 



- 

- 


- 

LS 5039 

0202950201 

Q 

1.14E-01 

- 

- 

< 4.12E - 02 

< 0.08 

- 

- 

- 



- 

- 


- 

LS 5039 

0202950301 

Q 

8.41E-02 

- 

- 

< 8.39L - 02 

< 0.28 

- 

- 

- 



- 

- 


- 

4U 2206+54 

0650640101 

Q 

5.04E-02 

- 

- 

< 1.21L + 00 

< 0.03 

- 

- 

- 



- 

- 


- 

Cen X-3 
Cygnus X-1 

0111010101 

0202401101 

I/E 

Q 

1.24E-02 

1.54E-01 

6 41+0.01 
-Q.Qi 

6 59+0 05 

42 06^'’-^° 
0.00 

5.77£ + 01!“2[ 
6.22£ + 00!™“ 

0 14+0.02 

0 01+0.01 

16.08 

3.46 

6 42+0 06 

6.50+°°^ 

qo + 100.30 

1.66L 

+ 

n7+8.33£+01 

^^-1.26£:+Q2 
a9 + 1.19£:+02 
^^-9.09£+01 
a9 + 1.65£:+02 
^^-1.18£+Q2 
a9+2.45£:+02 
^^-1.67£+02 

0.44+° 

0.46+° 

09 

18.93 

Cygnus X-1 

0202401201 

Q 

1.39E-01 

- 

- 

- 

- 

- 

944.58+'f,^5 

1101.52+lfl 
1243.17+1™ 
1199 67+^+™ 

i- lyy.Kj I 

1297.32+®+™ 

1.31L 

+ 

19.49 

Cygnus X-1 
Cygnus X-1 
Cygnus X-1 
Cygnus X-1 

0202760201 

0202760301 

0202760401 

0202760501 

Q 

Q 

Q 

Q 

2.94E-03 

5.82E-01 

6.18E-05 

1.38E-01 

6.54!;| 

6.48!||^ 

^■'^'^-0.02 

0.00 

0.00 

0.00 

0.00 

1.07£ + 01!J™ 
5.54L + 00 !™ 
6.01£ + 00!3;S 

1 1 9 F j_ m +5.53£'4-00 

0 01+'’ ™ 

0 ni+o+^o 

0 01+0.00 

0 01+0.00 

5.82 

7.72 

7.47 

8.97 

6 41+5.04 
^■^^-0.04 

6 47+0.03 
-0.03 

6 45+5.03 

u.j)u_0 22 

4.82L 

3.29E 

4.64E 

6.03E 

+ 

+ 

+ 

+ 

0.65™ 

0.61! 

0.61! 

0.83™ 

0.09!° 

04 

05 

04 

06 

02 

11 

03 

5.82 

> 25.5 
4.02 

> 25.5 

Cygnus X-1 

0500880201 

D 

6.15E-02 

- 

- 

- 

- 

- 


2.13E 

+ 

Ai +1.59£'+01 
^^-1.10£+01 

4.99 

Cygnus X-1 

0500880201 

Q 

3.73E-04 

6.41 

0.00 

4 Q4F + (-)(-)+ 1.25£:+00 

7.05L + 00_2 9g£^oo 

0.02 

6.95 

- 

- 



- 

- 


- 

Cygnus X-1 

0610000401 

Q 

1.43E-01 

^■^^-0.09 

82.681^0 

0 

^■^"^-0.01 

5.09 

- 

- 



- 

- 


- 

AXJ1749.1-2733 

0510010401 

Q 

6.80E-01 

- 

- 

< 6 .11L + 00 

< 0.12 

- 

- 

- 



- 

- 


- 


Table A. 3; Fit parameters of FeKo- and broad Fe feature detected. We present the centroid energy (keV), the width (eV), the luminosity (photons/s), EW (keV) and the statistical 
significance of the Gaussian lines. In the cases where we don’t detect FeKcr (cTsign < 2), we give the upper limits of EW and the luminosity. Parameters frozen or unbounded 
are included without an error estimation. Therefore they are not used in the plots and the subsequent discussion. The possible states are: quiescence(Q), flare(E), eclipse 
ingress/egress(I/E), eclipse(E) and dip(D). 












Appendix B: Spectral atlas 
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Fig. B.l: BeXBs data, model, model components and ratio data/model. The spectra are typically soft, with no Fe emission lines. 





































































































































































































SGXB 



Fig. B.2: continued. 







































































































































































































































































































































































































Fig. B.2: SGXBs data, model, model components and ratio data/model. The spectra are characteristically affected by high absorp¬ 
tion, with Fe fluorescent lines. 









































































































































































































































































































































































































































SFXT 












Fig. B.3: SFXTs data, model, model components and ratio data/model. 







































































































































































































































































































































































































7 Cassiopeae-like 



Fig. B.4: 7 Cass-like systems data, model, model components and ratio data/model. The data is usually fitted using thermal models, 
including Fe recombination lines. FeKo- is also usually visible. 


























































































































































































































































































































































































































HMGB 



Fig. B.5; HMGBs data, model, model components and ratio data/model. Soft spectra with no sign of Fe lines. 


Non Classified 



Brmiff, (taV) 


Fig. B.6: AX J1749.1-2733 data, model, model components and ratio data/model. We have no references for the luminosity class of 
the optical start, although the high absorption observed points to a supergiant companion. 





































































































Peculiars 



Fig. B.7: Peculiar sources data, model, model components and ratio data/model. These sources can be hardly categorized in any of 
the described HMXBs standard groups. 












































